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Abstract 
The Mississippian Exshaw Formation of Alberta consists 
largely of black, and in places calcareous, radioactive 
Shales. Forty-two core samples of the formation have been 
analyzed for some forty elements to quantify the naturally 
occurring radiosotopes U, Th and K and to study the 
geochemistry of these uraniferous black shales. 

Uranium concentrations up to 92.2 ppm (twenty-five 
times the average shale value) are largely responsible for 
the anomalous radioactivity, though in places potassium 
concentrations up to three times the average black shale 
concentration, are also a major contributor. The median 
Exshaw Shale thorium concentration of 9:2 + 1.9 ppm is in 
agreement with the relatively constant average shale Th 
values. 

Uranium is found to be correlated with organic carbon, 
suggesting that following transportation as a soluble (most 
likely carbonate) complex the uranium was precipitated and 
adsorbed by the carbon in a reducing environment. 

Thorium shows strong correlations with most of the 
hydrolysate elements suggesting that the element is bound 
mainly to clays. 

Compared to average shales the median Exshaw Shale is 
enriched an U, V, Zn, Sb; Pb; Ni, Mo, Cu and As while Cr, 
Ga, Mn and the rare earth elements are depleted. The Exshaw 
Shale rare earth element patterns are very similar to those 


of the North American Shale Composite and the post-Archaean 
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Australian Shale except for the presence of a negative 
cerium anomaly. This anomaly is thought to be due to 
dilution of the black shales by calcareous material and a 
Siliceous biogenetic component. 

A number of instrumental methods, including 
instrumental and epithermal neutron activation analysis, 
delayed neutron counting, fission track analysis, X-ray 
fluorescence and gamma-ray spectrometry, utilized for 
uranium, thorium and potassium analysis are compared. The 
relative merits and drawbacks of each technique are 
presented in an attempt to aid the potential researcher in 
deciding which method is most Suitable for his needs. A 
Vargationeofmthe “traditional method of uranium. thorium 
and potasSium analysis via natural gamma-ray spectrometry is 
presented and shown to require less sample and to be more 
rapid than conventional gamma-ray spectrometry, but with 
comparable accuracy and precision being obtained. 

The accuracy, precision and multielement nature of 
instrumental neutron activation analysis 1S exemplified by 
the analysis of a number of the more traditional USGS rock 
Standards together with three of the 'new' standards, SCo-1, 
SGR-1 and MAG-1 (two shales and a marine mud respectively). 


The results for these new standards are presented. 
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Tew introduction 
Since the 1890's it has been recognized that marine black 
shales often contain anomalous concentrations of uranium 
(Bell, 1978). tn Alberta a number of horizons are 
anomalously radioactive and are used by the oil and gas 
industry as marker horizons when logging and correlating 
wells. Examples of such units within Alberta include the 
Exshaw Formation and the Nordegg shale of the Fernie Group. 
Until recently however, few data existed on the uranium or 
potassium contents of the Exshaw Shale (Campbell, 1980; 
Havard, 1967) and to date, no data have been published on the 
thorium contents of the shale. 

The initial objective of this thesis was therefore to 
quantify, and if possible explain, the uranium, thorium and 
potassium contents of the anomalously radioactive pelites of 
the Exshaw Formation of Alberta. As this in part required 
the setting up, and in certain instances development, of a 
variety of methods for the analysis of U, Th and K (for 
example natural gamma-ray spectrometry, delayed neutron 
counting and fission track analysis) it became apparent that 
a description and comparison of the various techniques 
utilized would be both a useful and valuable contribution. 

Forty-two core samples of Exshaw shale, taken from 19 
wells, were analyzed largely by instrumental neutron 
activation analysis (INAA) and X-ray fluorescence (XRF), for 
some thirty-five elements in total. The locations of those 


wells sampled are given in Appendix 2. Samples were selected 
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from core areas of maximum and minimum total radioactivity 


using an Urtec minispec with a small Nai(Tl) crystal. 
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2. Neutron Activation Analysis 


2 eI Nt roducti on 

The principles of neutron activation analysis (NAA) 
were first developed by Von Hevesy and Levi(1936). However, 
1t was not until the development of the nuclear reactor that 
the technique could be used in a routine manner. Over the 
last four decades the increased availability of 
research-type nuclear reactors, together with the 
development of semiconductor detectors and the complementary 
solid-state electronics capable of collecting and reducing 
complex spectra have made NAA an extremely attractive 
method. At present the sensitivity, precision, generally 
non-destructive and multi-element nature of NAA makes it a 
competitive means of analysis when access to a nuclear 
reactor is available. A brief introduction on the theory of 
NAA will be given so that the reader might appreciate some 
of the later discussions. For a more complete coverage of 
the technique and its applications see Rakovic (1970), 
Keugere: (971)) DerSoeteset al u( 1972); aMueekesetyaie (1980), 


and Amiel et al. (1981). 


2.2 The theoretical background 

Unlike many techniques commonly used in the earth 
Sciences, NAA is based upon modification of the nucleus, not 
the orbiting electrons of the atom. Nuclei of stable 


isotopes are irradiated with neutrons which, via neutron 
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Capture, produce radioisotopes. These radioisotopes then 
decay according to a well defined time constant (termed the 
half-life, T1/2). The decay of many such radionuclides is 
accompanied by the emission of gamma-rays characteristic of 
the radionuclide. Detection and counting these gamma photons 
affords a means of elemental analysis, the energy of the 
gamma peaks being used to identify the radioisotope. The 
photopeak intensity is related to the element mass by the 
activation analysis equation:- 


MeN. OG. 0.1. 
A =. (1-e “Ati).e-Atd.(1-e-Atc) 
M.A 


photopeak area 

mass of element 

Avogadros number 

fractional isotopic abundance of isotope 
molecular mass of target element 

decay constant (0.693/T1/2) 

effective cross section 

HeveRoOnet Um Vnechae Sey) 

fractional detector efficiency 
fractional gamma yield 


HonoeoOoreztozs YP 


Ht th Uy TTP a 


tl lrradiation time in seconds 
td decay time in seconds 
ee count time in seconds 


Within a reactor there is a spectrum of neutron 
energies. They range from the highly energetic fission or 
fast neutrons produced by the fission process, through 
resonance or epithermal neutrons (also called intermediate 
neutrons) which are in the process of 'slowing down" due to 
elastic collisions with the reactor moderator; down to the 
lowest energy thermal neutrons. It is the thermal neutrons 
which are largely responsible for the (n,y) reactions. Fast 


neutrons, due to their greater energy, may induce reactions 
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which result in the emission of a charged particle (such as 
a proton or a-particle) or another neutron from the excited 
nucleus. The (n,p) and (n,a) type reactions are known as 
transmutat ions as the atomic number of the product nucleus 
differs from its stable parent. Although analytically useful 
in a number of instances, such reactions are often the 
source of interferences (see page 28). 

The neutron flux ¢, is the product of the neutron 
density (n.cm-°) and their average velocity (cm.s~') hence 


Givingtathesunits neemas is Due to the different average 
velocities and densities of the thermal, epithermal and fast 
neutrons there is a particular thermal, epithermal and fast 
neutron flux for any particular reactor operating level and 
position relative to the core (the source of the fast 
neutrons). Flux monitors of Au, Mn and Ni for example, may 
be used to meaSure the various components of the total 
neutron fluxes (De Soeteset aliigp.70-72, 1972))s. 

The probability of a particular nuclear reaction 
occurring is known as the reaction cross-section, and is 
given the symbol, 6. This is an energy dependent function 
ande1ssmeasuredain, barns, (iebarn =)1044 cm?) ssPractically 
all fast neutron reactions have very small cross-sections in 
the order of millibarns (mb) while the majority of thermal 
(n,y) reactions range from near unity to several thousand 
barns. The effective or reactor cross-section for many (n, 
y) teactions includes a component due to the resonance or 


epithermal neutron flux. This.effective cross-section, for 
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any particular (n,y) reaction, is given by:- 
gepl 
6.CLLeCtive = 68. 1) Les. 
gth 
Wheneriresy = 31... °2 0145-6 

The resonance integral I, and thermal neutron 
cross-section 6,, are available in tables (e.g. Lederer et 
al. 1978, Gryntakis and Kim 1976, 1978) while the epithermal 
to thermal neutron ratio is dependant upon the individual 
reactor and irradiation position and, as noted earlier, may 
be measured via flux monitors. 

The activity of a radioisotope following a ti second 
irradiation is governed in part by the isotopes half-life 
and 1S given by the growth factor 1-e°-Ati. The decay factor 
e-Atd, 1S a measure of how much of the initial activity 
remains after a decay time td. If after this time the sample 
GS COUNTeEd foresayperiod tc, the fraction.of the integrated 
area under the decay curve is given by i-e-Atc/ (assuming 
all emissions are detected). 

The number of gamma-photons emitted per disintegration 
of a particular nuclear species is known as its fractional 
yield or y-branching ratio and 1s given. the symbol 1. 

The area of the photopeak however, depends upon the 
fractional efficiency of the detector, e«. This is controlled 
by the sample geometry (i.e. the solid angle subtended by 
the detector) and the energy of the emitted radiation and 
fanaliy upon the intrinsic ett icrency or thesderecror 


itself. Efficiency curves for the various counting 
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geometries are uSually found empirically from the analysis 
of mixed radionuclide reference sources (e.g. ‘*?Bu, **Y, 


SOL ahdt Sb) Motaknownwacrivacy, 
2.3 SLOWPOKE II Nuclear Reactor 


2, 3,408 Introduction 

A supply of neutrons is essential for delayed neutron 
counting, fission track and neutron activation analysis. 
Invariably, a nuclear reactor is used to generate the high 
neutron flux necessary for practical analyses. 

As all three of the above methods of analysis were 
Carried out using the SLOWPOKE II nuclear reactor at the 
University of Alberta, it is felt pertinent to include in 
this thesis a description of the SLOWPOKE and some of its 
more important features. Most of the following information 


al.(1973) and the Atomic Energy 


is synthesised from Kay et 
DenCaonadashimicted WARCL) Iaterature se the reader ase retered 


to these publications for additiional information. 


2.3.2 Fundamentals of the SLOWPOKE II Nuclear Reactor 

The reactor was developed by AECL at the beginning of 
the 1970's to be used at universities and hospitals for 
isotope production and neutron activation analysis with a 
minimum of supervision. The word SLOWPOKE is an acronym for 
Safe LOW POwer Kritical Experiment. A maximum thermal 


nmegerron tlux ote» 10 “nm cms - US attainable, which as 
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greater than can be supplied by most neutron generators or 
sources, (forvexample = 94¢f):,, but at fer fessecost and 
Supervision than conventional reactors. 

SLOWPOKE is a small pool-type reactor (see figure 2.1) 
With a critical mass of enriched (to 93%) #7°U of 
approximately 0.85 kg. The small critical mass is the result 
of the beryllium annulus, bottom plate and upper shims which 
surround the fuel cage (see figure 2.2). Due to its low 
atomic mass, Be reflects a certain percentage of neutrons 
back into the reactor core where, once thermalized, they may 
bartictpacesrurther on themfission of 24'tL, 

The insertion or removal of a central cadmium control 
rod 1s used to control the fission process. This is due to 
the very large thermal-neutron capture cross-section of 
im Cas(cas 207000, barn). 

In addition to acting as a coolant and moderator, 
normal (i.e. light) water is one of the primary safety 
features of the reactor. Due to a loss in density as the 
water temperature rises its ability to moderate the fission 
neutrons decreases. The water is said to have a 'negative 
temperature coefficient of reactivity'. Consequently, fewer 
fission neutrons are thermalized and the chain reaction 
slows and may reach the point were it is no longer 
Self-sustaining and becomes Sub-critical as the fission of 
2331] occurs dominantly with thermal neutrons. Consequently, 
if there was a malfunction in the water cooling system the 


reactor power level would fall to where the amount of heat 
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being generated by the fission process equalled that being 


lost touthe= pool) surroundings: 


225oa)rocadiation: facilities 

The University of Alberta SLOWPOKE reactor has five 
inner irradiation sites (within the Be annulus) and one 
outer irradiation site (outside the Be annulus); see figure 


2,2scamples are =transtered rapidly (ca.) 1.5°s* one-way) 
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Figure 2.2 SLOWPOKE reactor critical assembly and core 


to-and-from the sites pneumatically via a 'rabbit' system. 
The inner site vials are 7 cm*® in volume and may accomodate 


two smaller 1.5 cm* vials while the large outer site vials 
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The U-Al alloy fuel rods are arranged so as to produce 
an equal, homogeneous neutron flux at each of the five inner 
Sirece Ryangeread ss io7 oe) mcarmied out a number of flux 
measurements on the Dalhousie SLOWPOKE reactor. From 
measurements on three inner sites the neutron flux was found 
to differ by < +2% (two sigma). Further measurements showed 
that the flux homogeneity varied <1% cm™' vertically in both 
the inner and outer sites. Radial homogeneity was found to 
Vary <1% em “in the inner site but waried by 5% cm * wn the 
outer Site. From their research they concluded that the 
neutron Hbux:= 

"..1S Stable, homogeneous and reproducible 

from day to day over a period of months." 

Ryaneeieale tpg 7s (OV Ge 
Similiar stability was found by Bergerioux et al. (1979). 
This is extremely important as it makes the simultaneous 
irradiation of flux monitors and/or standards with samples 
unnecessary. 

At the University of Alberta SLOWPOKE reactor the 
thermal to epithermal neutron ratio was measured via Au and 
Mn flux monitors (De Soéte p.70, 1972) and found to be 16.9 
+ 0.8 for the inner sites. The thermal to fast neutron ratio 
fOuRUNCmInne GEST tes «is 6400 a nO 922 a Ryan tall nor 
Zikovsky and Galiniver, 1984). As sis| Shown later the 


epithermal and fast neutron components may be used 


advantageously for the analysis of certain elements. 
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2.4 Instrumental Neutron Activation Analysis of Samples 
Two irradiation schemes were developed to maximize the 
number of elements determined together with a reasonable 
sample throughput. They included:- 
a) 4 Minute arradiation ati2.5 x 10 @n om 1s 
Db) 2 snovunrieradsavioneaty lo Oexe 0c naem- sr 
Counting was Carried oll with am Ortec 80 cm? active 
volume Win'® coaxial Ge(Li) detector, coupled to a Nuclear 
Data (ND) 660 multichannel analyzer with an Ortec 572 
amplifer and ND 575 ADC. The detector specifications include 
a relative efficiency of 18.5% and a measured FWHM of 2.1 
keV, sand peak=to=Comptonr ratio of 5341 for the 1332 kev 
photopeak of *°Co. Signals were assigned to one of 4096 
channels and the spectra collected stored on floppy discs 


for later analysis. 


2.4.1 Short Irradiation (4 minute) 

A number of elements, including Al, Mn, Mg, Ti and V, 
only produce short-lived radioisotopes upon irradiation with 
thermal neutrons. For their measurement the samples must be 
counted soon after irradiation, before the activity from 
these short-lived isotopes decays away. Table 2.1 gives the 
relevant nuclear data for the short-lived radionuclides 
routinely determined in this study. 

Samples and standards (both liquid and reference 
Matemials)MWere Irradiated satwic, ok UG CN no ee COM nO Ur 


minutes in site 1 of the SLOWPOKE reactor. Each sample was 
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Table 2.1 Nuclear data for short-lived radionuclide analysis 
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counted twice, initially after a twelve minute decay period 
for 300s live time at 10 cm geometry and then for 600s live 
time at 1 cm geometry following a thirty minute decay 
period. This scheme allows for some 12-13 samples to be 
analyzed every four hours. 

Although the majority of the elements in table 2.1 may 
be determined in both counting intervals the first count 
produces more precise results for Al, Ti and V on account of 
their shorter half-lives compared to the longer-lived 
isotopes. After the thirty minute decay period, the ?*Al 
activity and its associated Compton continuum has decayed to 
approximately 1710000 of -Lts initial value. This results in 
am increase in the Signal—to-background ratio for the 
longer-lived isotopes which, together with the longer 
counting time, lowers the detection limit and increases the 
precision with which they may be determined. Figures 2.3 and 
2.4 show sample spectra abtained following the two decay 
periods for example with a number of photopeaks identified. 

Analysis was carried out via an adaptation of the 
comparator method of INAA (also called the 'semi-absolute' 
method by Bergerioux et al. 1979). Conventionally the method 
involves the simultaneous irradiation of a known mass of a 
particular element with a sample containing an unknown mass 
of the same element. If the counting geometries and 
trradtation-cecay-count histories are identical for both= rhe 
Samplesand standard, then thewactivities Of the 1sotope ot 


interest in the standard and sample are directly 
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Figure 2.3 The y-ray spectrum of a shale after a 12 minute 


decay period 
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Figure 2.4 The y-ray spectrum of a shale after a 30 minute 


decay period 
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proportional to the respective masses of the elements. An 
unknown mass may therefore be determined by ratioing the 
activities and multiplying by the known mass in the 
standard. As mentioned earlier, when discussing the salient 
features of the SLOWPOKE reactor, one of its assets is the 
Stability and reproducibility of the neutron flux over long 
periods of time. This makes the simultaneous or regular 
irradiation of standards with, or following samples, 
unnecessary (Bergerioux et al. 1979). 

Single element liquid standards covering a wide range 
of concentrations were prepared from a variety of atomic 
absorption standards (Aldrich, Alpha and Spex) for 
comparison. They were heat sealed in 1.5 ml polyethelene 
vials and analyzed utilizing the previously described 
scheme. From these runs the specific activities 
(counts/microgram for each nuclide) were determined for the 
two counting periods. To reduce the uncertainty in the 
specific activity determinations, the liquid standards were 
re-analyzed a number of times following suitable decay 
periods. This ultimately reduces the error when determining 
the elemental concentrations of the unknown samples (taking 
into consideration counting statistics of the individual 
peaks). The specific activities and mean detection limits 
for each element determined in the study samples are listed 
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Teablee2zv2 Spectiic activities and meanidetection limits for 


Element 


Short irradiation 


Net photopeak activity 


(cts/ug) 


Detection Limitt 
(in 100mg shale) 
ee ee ee ee 


Na 65.24 + 0.6? Cup es 9 

Al 2.654 + 1.0! 310 + 180 

Ca Ons67 tite 1900 + 1400 

K 1.56 + 0.5? 1550 + 1100 

Mg ORV 2 Doe 1300 + 800# 
Ti 1.88 + 2.5) 900 + 540 

Mn 5782 + 0.2? 0.44 + 0.30 

V 158.4 + 1.0! 105) 72.0 

U 17276 + 0.2? 0.4 + 0.3 

Sr 42.5 + 0.4! 90 + 57 

el 2.26 + 3.0? 50 + 30 

Co 1516 + 0.4? 4+ 1 

Ba 103.1 + 0.2? 45 + 27 

Cs 616.4 + 0.5? 7 +44 

Dy 30172 + 0.2? Onze Oat 

Sm 5980 + 0.2? 1.0 + 0.6 

Eu 2394 + 2.9? Ou 20 ne 

In 48954 + 0.1? 4.6x10-? + 3.0x10-° 
Si 0.0084 + 3.0? 2.4x105 + 1.2x10° 


—+ * 


2 


coefficient of variation in % 

mean detection limit in ppm calculated from 40 
shale samples (+ 1 std. deviation) 

from the 4:12: 5 min scheme at 10 cm geometry 
from the 4:30:10 min scheme at 1 cm geometry 


** ASSumming no ?’Al(n,p)*?’Mg interference 
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Dead time and Pulse Pile-up Effects 

Highly active samples (due largely to 7#Al, °*Mn and 
**Na for example) can introduce significant error when 
determining element concentrations due to dead time and 
pulse pile-up (also called ‘random summing') effects. As 
unknown samples and standards rarely, if ever, have the same 
activities, corrections must be applied. Consequently, a 
deaGdetimeéscorrections (after Takeuchipet ale 1980) for 
short-lived radioisotopes in the presence of a high, longer 
lived background (e.g. ?*Na and **Mn) was applied to each 
peak and a correction factor for random summing effects was 
calculated and applied following the procedure of Wyttenbach 
C1o9 19. 

Data reduction 

All analyses were carried out off-line, in batch mode 
using the peak search programme PREP10 (Apps, 1978). The 
programme calculates amongest other data the energy, net 
area, background, counts s~', FWHM and % uncertainty in the 
net area determination of each peak. The user may select a 
variety of parameters including the minimum size a 'peak' 
must be before being processed, the Gaussian sensitivity 
(i.e. how much a peak may deviate from an ideal Gaussian 
shape before analysis for multiplet or overlapping peaks is 
initiated) and finally the number of iterations or ‘peak fit 
attempts' that are applied to multiplets (see also ND660 
Operator's Instruction manual). Recently, the programme was 


made available on the University of Alberta computer (an 
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AMDAHL 470V/6). In the future this will greatly reduce 
analysis time. 

A programme waS written by the author in FORTRAN IV 
which given the sample mass, live and clock times, together 
with the identity, area and uncertainty associated with each 
peak (taken from the PREP10 output) corrects for dead-time 
and pulse pile-up losses and calculates the corrected 
counts, mass, concentration and 1 sigma standard deviation 
for each element measured. Where necessary the programme 
also converts major element concentrations to oxides. Where 
a particular radionuclide decays with the emission of two or 
more quantifiable gamma-rays (for example **Na, °**Mn and 
3®Cl1) the additional photopeaks have on occasion been used 
when calculating the element concentration. As each peak has 
its own associated uncertainty the best estimate of the mean 
value of the element and its associated uncertainty, was 
found by taking the weighted mean of the individual peak 
determinations. The reported mean values, determined from 


replicate analyses of USGS standard rocks, were calculated 


in the same manner. 


2.4.2 Long Irradiation (2 hour) 

Certain elements including As, Cr, Au, La, Fe and Tb 
only produce long-lived (T1/2 > iday) radioisotopes when 
irradiated with thermal neutrons. Other elements, such as 
Ce, Sc, Sb and Th, produce both long- and short-lived 


nuclides when irradiated. Invariably however, these elements 
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are determined via their long-lived products for convenience 
and increased sensitivity or precision. 

To determine a number of the longer-lived isotopes 
standards and the 100 mg samples used in the short 
irradiation were re-irradiated for 2 hours in one of the 
Bivevinner SLOWPOKE sites ata lfluxs of Gx 10M-necm 2 sy *, 
They were counted five times over the 21 days following 
irradiation in order to maximize the number of elements 
determined with optimium precision. The decay times, 
counting conditions and elements determined in each counting 
interval are given in table 2.3. The table shows that the 
concentrations of a number of elements are repeatedly 
determined via the same isotope. To reduce analysis time it 
would be possible to replace the first three counting 
intervals with a single one hour count after two days or 
just abandon the second day count altogether. Increased 
precision is however obtained by taking the weighted mean of 
these independant analyses. 

Sample throughput is limited by the 10 hour sample 
count approximately 3 weeks after irradiation. Therefore, so 
as not to monopolize the ND 660 spectrometry system this 10 
hour count was usually carried out overnight and only four 
samples were irradiated per week. Using approximately 0.75 
ml vials (cut-in-half 1.5 ml vials) 4 samples, each 100 mg 


mass, may easily be packed into one medium sized (7 ml) 


Irradiatvonevial. 
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Table 2.3 Elements determined and counting scheme employed 
LOluOwiInge2 nour lreradiat ion. 


Decay Counting Counting Elements 

interval goemetry* time Determined 
ia. 1O0eem: 20 min Na,K,La,Eu,Sm,As,Ga,Sb,U(Np) 
el 1-3 cm.  \eve Na,K,La,Bu,Sm,As,Ga,Sb,U(Np) 
See j=3eem: ihe uate Na;K,La,Sm,U(Np) ,As;Sb, Sc 
6 as i em. i igke SG, he) Costa Ce Gite e. 


U(Np) ,Th(Pa) 

Circ liGn tL rekaies TORhe SCe Rib oreo Bare se ConGi- Hie 
Ni(€o) ,Fe,Co,Se,Tb, Yb, hu, 
Ta,Th(Pa),Fe(Mn), + Nd 


* Initial 1 day count outside Pb cave, all other counts in 
(Oem Pb cave to reduce background. 


Data reduction 

As with the short-lived INAA, spectra were stored on 
floppy discs for off-line analysis. All spectra were reduced 
using the batch programme NAARUN (Apps and Apps, 1981). The 
programme is an extension of PREP10 (Apps, 1978) and is 
fully imteractive. In addition to the peak search output of 
PREP10, NAARUN identifies the isotopes possibly responsible 
for the emissions. Given the sample mass, counting time, 
plus the dates and times the samples were removed from the 
reactor and counting began, it calculates the mass (in ug) 
and concentration (in ppm) of each element identified. At 
present, the data reduction is done off-line on the 


mini-computer based ND 660 system. 
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It cannot be over-emphasised how important the 
development of this programme has been in making long-lived 
INAA at the SLOWPOKE reactor facility routine in nature. 
Using the programme analysis of a single spectrum takes 
approximately 10 minutes. Where a similar analysis attempted 
manually using the PREP10 output and pocket calculator it 
would require a minimum of 5 or 6 hours. Once on the 
University computer, analysis time will be cut to a fraction 


of the present ten minutes. 


2.5 Interferences 
Like most other methods of elemental analysis neutron 

activation analysis may suffer interferences which, if not 
eonsudered, «can cesult in significant error. Such 
interferences fall into one of two catergories, namely: 

1) spectral and 

2) mucléar. 
Table 2.4 lists examples of both interference types which 
may occur when analyzing geological materials via 


short-lived radionuclides. 


2.5.1 Spectral Interferences 

As the name suggests, Spectral interferences result 
when two or more photopeaks (produced by different nuclides) 
lie within the resolving power of a spectrometry system. 
Generally, with Gaussian curve fitting applied to spectra, 


the win’: Ge(Li) detector and Tts associated electronics can 


; - -_ 
a, 
- _ 7 
: ~ 
: , 83 Seetae Mi una. baetonzs ll 
| > iS as 
sestinencs po lfgon gee she Satie? art? cP a 
yyeai \t ys-vad WPT eget 2agsee% aca t 
oy vy t sis pen: 1) = bm PE. On py A paNe 
ine fu! : j ‘ave ts e0n i G! 7 
} ri: oe) piaain te 
. o aad met 7 - 
bd { ¢ re= 8 
4 
afren zis Tk 
i352 . cbiawicn 
- 
ANMeg ss Th 6 2 rd S , ~esebienas 
4 ¥ ae | > ee ee | aA J ‘ —_ * Ae | \ot ves 
i@ 
qe | on 
ug ty agge® wae rs oe Ls eugene. él 5.6.6 
7 Sy, ira vie) n@Aae wa a 
i aie Goud sq 
7 7 
a 
7 Th "yaa teal! teed 
a ee 
; _ a 
Teeth fines omy 5a" ni wag ow % —_ 


7 


pea taye ~4 
- : 


Tue ™ sa ; r 7) 


ZD 


Table 2.4 Potential nuclear and spectral interferences 


Potential Nuclear 
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Potential spectral 
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Satisfactorily resolve two gamma peaks whose centroids are 
Seno KEVGapart, **Sucheis thetcase with the ="Movand °°Mn 
photopeaks at 843.8 and 846.6 keV respectively following the 
12 minute decay period of the short irradiation scheme. 
Following the 30 minute decay period however, broadening of 
the 846.6 keV Mn photopeak (due to the increased photopeak 
area resulting from the longer counting period and increased 
counting efficiency) makes the determination of Mg from the 
843.8 keV peak difficult. In such instances, or where more 
than two peaks are involved, or two peaks are too close to 
be satisfactorily resolved, the analyst may use another 
gamma-emission for analysis (if available, of sufficient 
intensity and free from interferences). In the above Mn and 
Mg example the analyst may use the 1014.4 keV (28.5% 
intensity) and 1811.2 keV (27% intensity) gamma peaks of 
77™Mg and °*Mn respectively for analysis. 

Where the nuclide being interfered with decays with 
only a single gamma-ray (e.g. *'Cr and *?V) or with only one 
eMmisslomvol SUfETcrent intensity for quantification (e.g. 
is Balor *-K) 7 stwWO "Options are open to the analyst. Errstly, 
if the interfering nuclide decays with the emission of two 
Or more gamma-rays the ratio of the area of the interfering 
photopeak with the area of an interference-free 'monitor' 
peak may be determined by irradiating and counting a 


Standard of the interfering element and a correction factor 


‘This is in part a function of the photopeak energy as, with 
increasing energy the FWHM also increases and resolution 


deteriorates. 
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calculated. Regardless of the irradiation, decay and count 
history of this standard, the ratio of the photopeak areas 
will be constant. Consequently, the counts to be stripped 
from the ‘photopeak’ of the unresolved peaks may be 
calculated by multiplying the area of the interference-free 
monitor peak by the empirically determined correction 
factor. Alternatively, if the interfering isotope has a 
havi-lite which is "significantly shorter than that’ of “che 
nuclide of interest the analyst may, by extending the 
‘cooling’ period, reduce the interference to insignificance 
while still retaining a satisfactory level of precision when 
determining the longer-lived nuclide. Analysis of the rare 
earth elements (REE) may be used to exemplify this latter 
point. Due to the extremely similar geochemistries of the 
REE in practically all cases the presence of one REE implies 
the presence cf the other members of the group. Thus, the 
determination of Lu soon after irradiation via the 113.0 kev 
photopeak of '7’m Lu (T1/2=155d) is interfered with by the 
113.5 keV photopeak of ‘'7°Yb (T1/2=4.21d). Following a decay 
period of 30 days however, the activity of '’*Yb will have 
decayed. tom<t7mcttits initwal aceivity «whivlee nehGewa li be 
ApRsge SY Oe lEeS Tin le tedmactiv ity. ssimilanly, Grtetie 
interfering nuclide has a half-life much greater than that 
of an element of interest, the interference may be 
calculated and corrected for by counting the sample twice. 
From an initial measurement the unresolved peak area may be 


measured while from a second analysis (some time later when 
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the isotope of interest has decayed away to insignificance) 
the interference from the long-lived isotope may be measured 
and when corrected for decay, the interference may be 
stripped from the photopeak measured earlier giving the net 
peak area due to the short-lived isotope. 

Less obvious spectral interferences may also result. 
For example, where a gamma emission is greater than 1022 
keV, Single- and double-escape peaks may be the source of 
interferences. Once again however, the main photopeak may be 
used to monitor the effect if corrections are necessary. 

The coincident counting of two photons (gammas and/or 
Korays) ate hidh countrates, or that Gecay in cascade, will 
result in a pile-up or coincidence peak that may be an 
interference. In the case of the single photopeak at 320.0 
keV of °*'Cr, summing of the two '’7’m Lu photopeaks at 208.4 
and 113.0 keV would produce an unresolved peak at 321.4 keV. 
Such peaks are difficult to identify and correct for as the 
severity of the interference is largely related to the 
sample activity, matrix composition and various radioisotope 
decay schemes. The effect is therefore best avoided by 
keeping the spectrometer busy time <10% by varying the 


irradiation:decay histories and/or counting geometry. 


2.5.2 Nuclear Interferences 
Such interferences result when the radionuclide under 
investigation is produced via a reaction other than that 


under consideration. INAA using thermal neutrons largely 
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utilizes (n,y) reactions. Consequently (n,p) and (n,a) 
reactions are potential interferences. 

As a reactor flux is rarely completely thermalized, 
neutron activation by fast neutrons followed by the emission 
of a proton or alpha particle is in certain instances a 
favoured form of decay. As noted earlier when discussing the 
SEOWPORKRES reacvor@aliast ilux of approximately 2.2 x 10° 'n 
cite Ss 8) may be sobtained at’ full power. This 2s utilized 
favourably in some instances? however, in other cases 
corrections need to be made. 

Zikovsky and Galinier (1981) calculated interference 
factous (iP) forssomer50(n,p) and i00e(n ale reactions for 
the fluxeconedi tions occurring, rnsa SLOWPOKE 11 reactor. The 
interference factor (IF) may be defined as the mass (in ug 
for example) of an element which, via an (n,p) or (n,a) 
reactions produces a Speciviicractivity equal “to: that 
produced by an (n,y) reaction of 1 ug of the element of 
interest. A selection of transmutational reactions and their 
IF values of particular interest to the geologist are given 
in table 2.5. From these values it may be calculated that 
1000 ug) of Al would produce, via the reaction 77 Al(nyp) = 7Mg, 
the equivalent *’Mg activity produced by the irradiatzonvot 
185 ug of Mg via the reaction ?*Mg(n,y)*’Mg.* From the 
irradiation of Al standards of various concentrations when 


determining the specific activity for Al it was found that 


2in the determination of Ni via the reaction **Ni(n,p)**Co 
and Felyvia “‘Fe(n,p) “Mn for example: 
2 1000ug Al/5.4 = 185ug Mg (equivalent) 
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Table 2.5 Examples of interfering transmutations. 


Reaction IF (ppm) 
2*Si (n,p) 2%Al 690 
a Moe (np) tena 980 
2 ew tie) cee Mn Tox Oe 
pelea Wet eee MiG 5.4 
Bid Goegany ga ay I> See SHON: 
a) A? owe se) Sp 7 Gud 
Wier Xi enn) = A yssih 1600 
Seen eee Na 1200 
Connon). eoMn Soe an KO) 
SoS ei Crew) a Ma 6000 
= Mine red ey door ke) 
SOC ren oe Alex 20s 
12?Gd (no) 72 Sm 4.9 x 105 


(taken from Zikovsky and Galinier, 1981). 


for the 12 minute decay period at i0 cms geometry: 
77™7Mg cts at 843.8 keV photopeak 
AANA = 0.16 + 7.2% (N=7) C1) 
MORAL 
and for the 30 minute decay period at 1 cm geometry. 
77Mg cts at 843.8 keV photopeak 
$$. = 0.71 + 0.8% (NAS) (2) 
ug Al 


while the specific activities for the Mg (n,y) reaction 


following the short irradiation scheme are: 


0.761 cts / ug Mg + 2.1% (N=5) at 10 cms geometry (3) 

4.04 cts / ug Mg + 0.8% (N=5) at | cm geometry (4) 
Thus from equations (1) and (3) 1000ug Al would give the 
eduivolent - Mg activity of 2 iGr2cug (7.5%) Mg while 
equations (2) and (4) give the equivalent of 176+4ug (2.3%) 
Mg. The weighted mean of these two results gives 181+4ug Mg. 
The agreement between the theoretical and empirical results 


is not quite as excellent as it first appears. Zikovsky and 
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Galinier (1981) used a gth/¢gfast ratio of 5 in their 
Galculations. In reality a value of 4.6 is more accurate. 
THis wouldschange the If for the 77Al(n,p)2@"*Morreaction from 
9.4 to 5.0 thus 1000ug Al would give a theoretical value of 
200 ug Mg. However, considering the theoretical nature of 
Zikovsky and Galinier's method, the agreement between the 
two results is still very good and may be used as a first 
approximation of the potential interferences. 

So that the geologist might more fully appreciate the 
possible extent of such interferences this author calculated 
the interfering effects of 1% of each of the ten major 
oxides routinely determined in geological materials from the 
IF values of Zikovsky and Galinier. From the results shown 
in table 2.6 one may conclude that the majority of the 
transmutations may be ignored unless rock-types of extreme 
composition are being analyzed. The 65-70% SiO, content of 
an average granitoid would produce the equivalent of c.525 
DomeAleCOS4% Al, Os) “whachtconstitutes=<i|yeoLrathestypical 
content of such rock-types. Conversely, the 40.5% S10, % 
content of the USGS rock standard DTS-1 (dunite) would 
produce the equivalent of c.320 ppm Al(0.06% Al,0;) which 
represents 25% of the recommended Al,0, content of the rock 
(02247 Al=O.)a0The Mg contentvotesuchsultramatic rock=types 
Mayeaiso result in Significant errors when determining Na. 
Still using DTS-1 as the example the recommended MgO content 


of 49.8% would produce the equivalent of c.300 ppm Na (0.04% 


* all recommended values taken from Flanagan(1976). 
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Table 2.6 Calculations of the interference effects of 
ba Of seach, Of the major cock forming oxides. 


Oxide (nos) (n,a) 

SauOy, 1.0 poem Al 0.9 ppm Mg 
Aly, 0; 980 ppm Mg 4,4 ppm Na 
T1075 O07 Sopn Sc = “eee 

Fe,0O, 0.04 ppm Mn <(i.i2e,ppm .Cr 
MgO 6.15 ppm Na 1.8 ppm Ne 
CaO CSmpia A eee 

Nese. © a. eee Sua DOM e 
K 20 0.24 ppm Ar ieee pom al 
MnO 645.5) DpmeCE 0.04 ppm V 
Peo. 4700 ppm $i 0.35 ppm Al 


Na,.Q) which is almost 600% greater than the recommended 
value of 0.007% Na,O. The MnO content of most rock-types is 
<1%, therefore the V interference is generally negligible, 
this is not always the case with regards to the (n,p) 
reaction. In the USGS standard W1 (diabase) the recommended 
MnO content of 0.17% would produce the approximate 
equivalent of 11 ppm Cr which represents almost 10% of the 
recommended 114 ppm. The importance of the Si production via 
2P.Criep) Si 1545 in cealdty,. vorelevant .as «the.Si icontent of 
samples when measured is determined via the reaction 

2 Si(oy oO) 2A at SLOWPORE.. 

Where necessary corrections for these transmutations 
may be made by monitoring the (n,y) photopeaks of the parent 
nuclides responsible.*® It is extremely important to note 
however, that if the analyst runs a Standard of a particular 
element to determine the extent of its transmutation 


5 see for example Kerr and Spyrou(1977) where the 1368.4 keV 
photopeak of 74Na is used to monitor and correct for the 
23Na(n,a)?°F interference to *°F produced via the reaction 
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imberterences the corrections obtained may only be used for 
eNesPanticular decay and) count histories that the 
corrections were determined for.* This is because the 
half-lives of the parent and the nuclide of interest are 
rarely if ever identical and the ratio of the monitor and 
peak of interest would be continually changing with time (as 


would be the correction factor). 


27 5.S uly) anterrerences 

The siaSsszoneproduets, of, 7°" U upon) which uranium 
analysis by delayed neutron counting and fission track 
analysis are based) are also a potential primary 
interference source. The fission of *°°U may be written as: 


735U + jn ——> 27340 —> xX +Y¥+ fin +9 


where: 
Q = the energy released per reaction (c. 200 MeV) 
f = the number of neutrons produced (c. 2.5/reaction) 
X = fission product nuclei 
Y = fission product nuclei. 


The two most probable masses of the fission nuclei are 95 
and 139, about the maxima of the bimodally distributed 
fission yield curve (figure 2.5). Following the method 
outlined in De Soete et al. p.481 (1972) correction factors 
fOr thestission of °° "°U in natural uranium were calculated. 
Those that were significant when analyzing geological 
materials are listed in table 2.7. The fission yields, 


thermal and resonance cross-sections, and remaining nuclear 


‘ unless the activities of both sample and standard are 
calculated to some reference time. 
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MASS NUMBER 


(Source: Katcoff, 1960) 


Figure 2.5 Fission yield curve of 7°5U induced by thermal 


neutrons 
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data were taken fron Katcorf (1960); Gryntakis and skKim(1976, 
1978) and Lederer et al.(1974). 

Most severely affected are the determinations of Zr, Mo 
and yBasvia the radioisotopes = °2r, @=Mo and #2°%Bayein 
addition, the light REE are also affected. This interference 
may be very important as it is the relative abundances of 
the REE that is of particular interest to the geologist. 
Consequently, care must be exercised when analyzing for the 
REE by NAA in the presence of uranium.’ If a uranium 
Standard is prepared and analyzed to measure the fission 
interferences empirically it is obviously very important to 
insure that the standard is of natural (or at least known) 
isotopic abundance. Empirical and theoretical (n,f) 
interferences were determined in this work. Their 
effectivness was tested whan analyzing the Canadian 
Radioactive Reference Ore DL-1 (U = 41 + 2 ppm) and the 
corrected results obtained for this standard were found to 
be in good agreement with those found by Jagam and Muecke 
(to be published, pers. comm. G. Muecke, 1982). 
Consequently, were necessary, these corrections have been 
applied to the samples of this study. While this thesis was 
in the latter stages of being produced a paper addressing 
the problem of 7*5U (n,f) interferences to REE analyses was 
published by Kennedy and Fowler (1983). Their findings are 
very similar to those measured and reported by this author 
Iie Le ei. 


7 100 ppm U in a sample would give the equivalent 
concentration of approximately 26 ppm Nd and 24 ppm Ce. 
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Table 2.7 Fission interferences due to the irradiation 
of natural uranium 


Isotope Fission Yield(%) lug U equivalent(ug) 
one cl Sow, 0.62 Ba 
peta 6.2 Gen Se Zc 
°°Mo 6.06 0.78 Mo 
ESA Wal == 0.02 La 
i Ce 6.0 0.24 Ce 
NG OARS) Oe26"Nd 
* produced during the decay of the fission product 
Pode vio Cer Teac tl Olle (Bo =———= 6 == elas lino = 12 fog) 
interference for decay time = 4 days. 


2.5.4 Other Interferences 

Secondary reactions induced by prompt gammas or charged 
particles emitted from the compound nucleus were considered 
but were concluded to have negligible effect when analyzing 
'typical' geological samples at the SLOWPOKE reactor. 

Similarly, second order reactions® need not be 
considered here as they uSually only present a problem when 
analyzing for a particular trace element in a near 
monoelemental sample that lies next to the trace element in 


the periodic table (see Kubota, 1977). 


2.6 Instrumental neutron activation analysis results 


2.6.1 Uranium, thorium and potassium 
Uranium analysis 
Invariably uranium has been determined via one or more 


Ofethe gamma-emissions of “-7Np as it decays tow] “Puscee 


"when the product of an (n,;7) ot transmutation reaction is 
itself further activated eg.*°Si(n,y)?'Si——> sp —~>’'P(n,y) 
ee peresultingin a.GainelieP Om Tal (n,y)vegad (ney ae la 
for a) lOSS OL an 1S50t0pe., 
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Mopgantand Lovering, 1963° “Schock, 1977* ‘Baedecker ettals, 
1977; and Hart et al. 1980). However, when ?°°U decays to 
272°Np it emits a 74.5 keV gamma photon of sufficient 
intensity to be used in quantifing thesuranium! comtent of 
many geological samples (see for example Vukotic, 1981). 

In this work uranium was determined from both the 
SHOntgrauee  ong-ivvede(n, 7) ®radvopnoduces ah 24 °u: 

The detection limit (Currie, 1968) £or uranium varied 
with the sample background when using the 74.5 keV emission 
of ?°°U. However, based on some forty shale samples, the 
mean detection limit for a 100 mg sample was 0.4 + 0.3 ppm U 
(at the 90% confidence level). 

Following thentwo hour wirradiation acraviilux of yl x 
10*29n em" “s-*sthe detection limit forsuranium, determined 
solely from the 106.1 keV emission of **°Np was found to 
talivtrom’ 2.3 220.7 *ppm Ustollowing acdecay perniodrot 
approximately e24 hours chow. 0e200.3 (ppmgUsaltern 2adays; «down 
tOLOnS ste0n2 ppmeUsatiter ea decay period of v6%days.SThis 
reduction is largely due to the simultaneous, but more rapid 
decay of 7*Na and its associated Compton continiuum. One 
major problem associated with the use of the 106.1 keV gamma 
emisSion is that it may suffer interference from the 103.2 
keV photopeak of ‘**Sm (T1/2 = 1.94d) even when a Ge(Li) 
detector and counting system of superior resolution is used 
(see Kuleff and Koslodinov, 1981) and Gaussian spectral 
EPtetinge 1s applied.” 


Use of a low energy photon spectrometer (LEPS) detector, 
with superior resolution to that of a Ge(Li) detector avoids 
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Similarly, the same interference to the 103.7 keV peak (K 

Re vays Of 7°" Pu) "makes this péak “offs lierle=suset ior analysis 
unless corrections are made. The 209.7 and 117.7 keV 
PHOLEOpeakSwOL <e=Nprabevinterferéea with bye ee Lu (91/2 155 
Guimande e Yb (T1/24=—30%6d) "sespectivelyasOr the 7277.5 and 
22626 KEV Dhotopeaks Of *°- Np the 277.5 keV emission 1S to 
be prefered as the 228.2 keV line suffers an interference 
prom the=229,3 «keV -photopeak of “8 2Taltht2'= 115d)ewhieh 
becomes increasingly significant the longer the decay 
peerod, sl neaddueion, «uhek-- “USrtissionwproduct “146 (hirssion 
yielde4.7%)ewith®a habf-lite, cfe7600 hours also emits ‘a 
gamma photon at 228.2 keV (G. Muecke, pers comm. 1982). 
Spectral interferences to the 277.5 keV photopeak were found 
to be negligible and the results obtained from it more 
accurate than those obtained from the other emissions. 
Unfortunately however, due to the combination of decreasing 
detector efficiency with increasing energy, and the lower 
intensity of this emission, the uncertainty associated with 
Eniselimewuseoreater ‘thanethat or ‘thes 06. for 22672 "key 
emissions. Under the present irradiation:decay:count scheme, 
andsusinguonly the 277 ,5ekeV emission, a concentration or 4 
ppm U in a 100 mg sample may be determined with a precision 
of +10 % or better (single analysis) following a decay 
period of 6 days. Conversely, the 74.7 keV photopeak of **°U 


may be regarded as interference-free unless samples 


*(cont'd)this interference but these detectors are generally 
not as efficient as the coaxial Ge(Li) detector and hence 
longer counting times are required for a precise analysis. 
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containing an appreciable amount of lead are to be analyzed. 
Upon excitation by more energetic gamma- and/or X-rays, lead 
may emit secondary fluorescent Ke X-rays with an energy of 
74.97 keV which cannot be resolved from the 74.7 keV gamma 


photons of 7°°U. Such samples were not found in this study. 


Thorium analysis 

Unfortunately, unlike 7**U, the sensitivaty forth 
analysis via the short-lived ?°*Th was found to be quite 
poor); (c. poG0lers/ugpThavsis J 7600s cessugel) se Fore thesshales 
analyzed, this resulted in a mean detection limit of 16 ppm 
Th andwaseduantatieatzon! dumitwot S0sppmethn (Currre, 1963). 
Consequently, although the 86.2 keV photopeak (2.7 % 
intensity) of 7°°Th was positively identified in a number of 
instances it could rarely be used for precise or routine Th 
analysis in this study. 

Following the two hour irradiation, and using only the 
311.9 keV gamma emission of ?*%Pa, the mean detection limit 
for thorium was caleulated? to be: 0.3 + 0.4 ppm followangi:a 
decay period of approximately 21 days. Based solely on 
counting statistics, 10 ppm Th in a 100 mg sample may be 
determined with an uncertainty of + 1% or better (1 sigma 
Standard deviation). As is shown later accurate and 
sufficiently precise results may be obtained following a 


decay period of 6 days. 


Potassium Analysis 
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Potassium was determined via its most intense gamma 
emission at 1524 keV (18% intensity). The emission is 
interference-free and because of the half-life of ‘4?K 
(12.2h) the same isotope is used following both the short 
and long irradiations. Following the short irradiation and a 
decay period of thirty minutes the mean detection limit for 
potassium (based on 42 shale samples) was calculated to be 
On 155% 22057105 COn1S6% +80. 132% K50) inca s00messamplel (90 
% confidence level). After the two hour irradiation the most 
precise and accurate results were obtained following a decay 
period of 3 days. 

Accuracy of the U, Th and K analyses 

To assess the accuracy of the U, Th and K analyses by 
both the short and long irradiation schemes, four 100 mg 
aliquots of the Canadian Radioactive Reference Ore DL-1 (U 
=41 +2ppm, Th = 83+5ppm and K = 2.50%) were prepared and 
analyzed in accordance with the previously described 
procedures. The weighted means, one sigma standard 
deviations, number of determinations and scatter of results 


obtained following the four minute irradiation are as 


follows: 
U = 40.8+0.1 ppm (N=4) scatter=6.54% 
Th = 79,.2+1.9 ppm (N=4) scatter=8.05% 
K = 2.48+0.07 % (N=4) scatter=2.50% 


Clearly the agreement between the recommended values and the 
results of this work are excellent. The scatter of the 


uranium and thorium results suggests however, that there may 
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be some sample inhomogeneity and/or that variations in the 
counting geometry have produced the scatter and may 
Enerefore antroduce a Significant source of random error, 

To further assess the four minute irradiation scheme 
for both accuracy and precision the USGS rock standards 
GSP-1, AGV-1, BCR-1 and W-1 were analyzed, the results are 
given in table 2.8 As thorium was only quantified in GSP-1 
values are not reported. As with DL-1i, the agreement between 
the recommended and measured concentrations are excellent. 

bollowing the two: hour irradiation of DL=1, and 
uti ligzingwonly: the277..9) 3112.9 and 1524 keV photopeaks of 
““eND,)- =Parand sack respectively, the individual results 
and weighted means of two separate analyses of DL-1 are 
reported in table 2.9 As can be seen the most accurate and 
precise results for potassium were obtained following a 
decay period of three days. Both the six day and ca. 
twenty-one day determinations of thorium are in excellent 
agreement with the recommended concentration. The uranium 
results obtained from the two 100 mg aliquots gave 
Signmiutacantly datferent results witha mean, Orwoo .020. 2 ppm 
(N=6) scatter=6.58%. As with the short irradiation results 
the scatter of the results from the long irradiation may be 
due to sample inhomogeneity and/or slight variations in 
counting geometry. Fission track analysis of an aliquot of 
DL-~1 would show whether inhomogeneity was a Serious problem. 

Finally, the standards GSP-1, AGV-1, BCR-1 and W-1) were 


also analyzed in keeping with the long irradiation scheme to 
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Table 2.8 Uranium and potassium results for selected USGS 


rock Standards following short irradiation 


URANIUM (ppm) POTASSIUM OXIDE ($%) 


5.50 + .07 (n=3) Sa) 
Aygo} &3 oik3} (le) 


1.88 (a) 
19 Bee Aan (T1=2) -10 (n=2) 2.89 (a) 
DOs 25 sks} (Ge) 


1.74 (a) 
1.69 + .24 (n=1) -16 (n=1) 1570) (a) 
1.81 + .09 (b) 


0.55 + .14 (n=2) 0.58 (a) -08 (n=3) 0.64 (a) 


(a) Flanagan (1976) 
(b) Millard (1976) 


Table 2.9 Uranium, thorium and potassium results for Di-1 


following a 2 hour irradiation and various decay times 
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evaluate the accuracy and precision of the U, Th and K 
determinations. The results, and a selection of recommended 
values, are given in table 2.10. Once again the agreement 
between the two is in general excellent. The uranium 
determinations are quite accurate though rather imprecise as 


a consequence of the low levels being analyzed. 


2.6.2 INAA of other elements 

Results of 4 minute irradiation of standard rocks 

tilizing the short irradiation scheme and data 

reduction method described the, USGS rock standards 
GSP-1(granodiorite), AGV-1(andesite), W-1(diabase), 
BCR-1 (basalt)a, PGC l(peradorite), DES71(dunate))) SCo—-1 Cody 
shale), SGR-1(Green River shale) and MAG-1(marine mud) were 
analyzeda The results obtained for GSP-1, DIS=1)7 SCo-1) 
SGR-1 and MAG-1 are given in tables 2.11 to 2.15 together 
with recommended and other literature values. As can be 
seen, the method is accurate and for most elements, precise. 
From the variety of geological materials analyzed the method 
is applicable over a wide compostional range. Its 
application to the analysis of ultramafic rocks is however 
limited due to the extremely low levels of many elements in 
Such rocks: 

Analyses of the new USGS rock standards SCo-1, SGR-1 
and MAG-1 show that acceptable results may be obtained for 
argillaceous materials and, in addition, adds to the limited 


data base on these materials (compared to the ‘traditional’ 
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Table 2.11 Short INAA results for USGS Cody Shale (SCo-1) 


Element This work Literature 
/Oxide (+ 1 std.dev.) valuest 
Al203 OSes oer Ole? US rat 5) 8 oesul hee ecrrel eS gel Ss 


MgO joo! TOhehe SrsCLURUO enoUL0 24, 2.00 

CaO 2 Oise te 08/0 Zhe OF 0,URle 2h n0Se 

Na20 OnS2e 25 07.04 ORe SD Ors 7) Ol Sle 0 O00 ano Oia O77 
K20 Cin] On tenia lke 2.6540.04',2.67540.0177,2.7°, 1.795" 
TiO? 0.54 + 0.04 ORG 2OROet = (0) 64-25 00.54 4.2 

Ba DiGoscts 2 G22 ay 4e6 226 560* 544754598" 

Gi 46 + 9 Sheu) 

Co Gi 0e te OLY Slee es Coes) 

Cs Set le eS bes) 6. 82207168 7/6". :G0 - 

Dy 34/7 te ONE ARAM eNGr toh Oe 7S 6 OOM 7 ae)” 

Eu doaese (te 0790-4 Teee2 Ee ORNs oy lee Oe EO 02 - aledon ea lee a. 
in Os Oot cee Cx Oe | Cee Pian se are 

Mn 36 Oia e.6 SIO tore Cole 

Sm Dewalt Ondo Dee SLO ears O40 5.07 = 

Sr poe Se 19 HOSe or pba se 2, 200s 17 

U S50 te UR 07 Se Ot OU Oem Oo 

V 1o2e te 3 PoGeeer si 20) a ciston 


(Major oxides in weight percent, trace elements in ppm.) 


+ References. 


' Fabbi and Espos (1976) 2 Thomas et al. (1976) 

> Katz and Grossman (1976) ‘* McLennan and Taylor (1980) 
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(references continued on table 2.11) 
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Table 2.12 Short INAA results for USGS Green River Shale 


(SGR-1) 
Element This awork Literature 
/Oxide Ce ed stg sdey..) values 


NLROx SATU O HS PE yaa at 
MgO 4.05 + 0.14 4.504+0.35',4.3240.027,4.5° 
Cao NS?) es Olan ene Chaco hon 
Na20 3.19 40.08 “2.6620.06",3.0340.017,3.0" 
K20 (66) 00.05  “f- 7 120.050) 1.6040 -012,1.6° 
TiO, 0.40 + 0.03 0.36+0.01',0.26° 
Ba 308 + 12 328+8' ,225+120°,278! 
Cl 39 + 27 444) 
Co Ged te a0 10.3+0.57,12.05 
Cs a awe Bs 4.6140.36°,5.16! 
Dy aes OMe Ti nee 
Bu 0.52 + 0.04 0.52+0.07°,0.50*,0.505 
In neas 0 ae aa sa 
Mn 250 + 8 251+12,232° 
Sm PUR (GL Piles einer: 
Sr 397 + 24 445+6',332+10' 
U 5.66 + 0.09 5.60+0.17°,4.85* 
Vv 130 + 3 13443" 


(Major oxides in weight percent, trace elements in ppm.) 
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Table 2.13 Short INAA results for USGS marine mud (MAG-1) 


Element 
/Oxide 
A120; 
MgO 
CaO 
Na;0 
K20 
T1i02 
Ba 
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Co 

cs 

Dy 


Eu 


Sm 


Spe 


This work Literature 
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(Major oxides and Cl in weight percent, 
trace elements in ppm.) 
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Table 2.14 Short INAA results for USGS Durnite (DTS-1) 


Element This work Literature 
/Oxide (+ 1 std.dev.) values 
A120; U2 292) 07027, Caras 

MgO 4875 2 0.9 aO".8 
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(Major oxides in weight percent, trace elements in ppm.) 
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Table 2.15 Shore uNAA eresults for USGS Granodiorite -(GSP-i) 


Element Thasiwork Literature 


7VOxide C2 ist cardev.) Values 
See ee ee be ee ee ed ed Oe a ee 
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Standards). 

rhe four Minute: irradiation resuits obtained for the 
Exshaw Shale samples are tabulated in Appendix 1 together 
with the results of the long irradiation scheme. 

Results of 2 hour irradiation of standards rocks 

Following the 2 hour irradiation at full power, element 
concentrations were determined via a variant of the 
fabsolute method’ of activation analysis. This method 
normally requires the determination of the neutron flux and 
Pes TOUSEriDULLOn vie ea flux moniltor(s), ipradiatred 
Simultaneously with samples. This, as is apparent from the 
earlier discussions on the SLOWPOKE II reactor, is 
unnecessary aS it has been shown by various workers that the 
neutron flux of the SLOWPOKE reactor is both stable and 
reproducible over extended periods of time. With the 
knowledge of the neutron flux, the mass (and hence 
concentration) of an unknown element may be determined from 
the activation analysis equation (page 5). Although in most 
cases precise, the method is often inaccurate due to 
uncertainties in the nuclear and isotopic data (e.g. 
half-life. cross section etc.). Analysis of the USGS rock 
Standards GSP-1, AGY-1, BCR-1, PCC-1 and DYTS-1 showed this 
to be the case for a number of elements. Consequently, in a 
manner similiar to that employed by Potts et al. (1981) 
correction factors, based on the average discrepancy between 
the 'found' and recommended values for the standards, were 


calculated. To test the method with corrections, samples of 
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the old USGS rock Standard W-1, and two of the new standards 
~CO-1 and SGR-1l)) were analyzed.*® The results obtained are 
bepoOGred itatables, 2.16: to 2.18 and as with the other 
Standards, recommended and published values are given for 
comparison. In practically every case there is very good 
agreement between the results of this work and those of 
other workers. 

Using these correction factors therefore, the Exshaw 
Shale results following the long irradiation, were 
corrected. These results are tabulated in Appendix 1 


together with the short irradiation results. 


2.7 EPITHERMAL NAA 


2 cio Incroduct2on 

AS mentioned in an earlier section, the effective or 
reactor cross-section for any particular (n,y) reaction is 
made up of two components. These being the thermal neutron 
cross-section 6, and the resonance integral I, (see page 7). 
With increasing neutron energy the majority of (n,y) 
reactions follow the so-called '1/V law', that is, the 
thermal cross-section is inversely proportional to the 
velocity (and therefore energy) of the incident neutrons. 
Examples of such reactions include ?’Al(n,y)7°*Al, 
23Na(n,y)?4Na and **Mn(n,y)*‘Mn. Certain isotopes however, 
with large 1,/G ratios, are preferentially activated by the 


‘'°These samples were not included when calculating the 
“eorrection.. factors. 
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Table 2.16 Results following 2 hour irradiation of USGS Cody 


Shale (SCo-1) 
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fable 2.17 Results following 2 hour trradiation of USGS 


Green River Shale (SGR-1) 


Literature 
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Table 2°16 Results following 2 hour irradiation, of USGS 


Diabase (W1) 


Element This work Literature 

/Oxide (= *T"std.dev.} values 
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(Major oxides in weight percent, trace elements in ppm.) 
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higher energy neutrons in the epithermal or resonance 
neutron range. By ‘filtering out’ the thermal neutrons; so 
that only the epithermal and fast neutrons may activate the 
sample, one may preferentially enhance the 
Signal-to-background ratio of those isotopes that have large 
resonance integrals over those which follow the 1/V law. 
This method or technique is known epithermal neutron 
activation analysis (ENAA) and may be used advantageously in 
certain instances as it:- 
1). lowers the detection limits for those elements which 
have isotopes with major resonance peaks (e.g. U and Th, see 
Gladney et al. (1978) and certain REE, Baedecker et al. 
(1977)) by reducing the background due to the major Compton 
producing isotopeswe.guc“Al,s2 ‘Na, > °Mneand: 4*Se) (common 
constituents of many geological materials). 
2). reduces the sample activity which: 

a) makes post-irradiation handling safer. 

b) results in shorter decay times being necessary 

following irradiation. 

c) allows for more sample to be irradiated if 

necessary (which allows a more representative 

Sample to be taken and thus lowers the detection 

limit for many elements). 

a) reduces corrections resulting from highly 

active samples (e.g. pulse pile-up and dead-time 

eorrectionsy. 


3). reduces the 7°°U (n,f) interference to certain isotopes 
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6.G. — Mo, "**Ba and some of the REE. 


2.7.2 Previous work 

Due to their very large thermal neutron absorption 
cross-sections, cadmium and boron have been utilized by a 
number of workers in ENAA (see Steinnes (1970), Bem and Ryan 
(igei), Bhmann et ale (1980) “and Gtuart and Ryan (1931))- 
MneweEOSS Cress -sections fore = Cdvand - 2 plouredsagainse 
meutron energy are shown in figure 2.6. At approximately 0.2 
eV’ Cd hasan absorption cross-section of about 6,000 b. 
Above thisjit falls off rapidly to =10 b between 1 and: 10 
eV. This abrupt change is known as the 'cadmium cut-off f'. 
Boron-10 on the other hand rigidly follows the 1/V law. At 
first sight cadmium appears an ideal thermal neutron filter 
for ENAA. Unfortunately however, the metal is toxic, thus 
fabrication of shields must be carried out with care. in 
addition, upon irradiation Cd produces a number of 
radioisotopes. One of the radioproducts, ‘‘'°Cd (T1/2=53.5 
hrs), has a high specific activity which makes handling of a 
shield soon after a long irradiation undesirable. 
Conversely, *“°B activates via the reaction *°B(n;a) 7bi 
producing stable nuclei (see under DNC) resulting in a 
Minimum of Shield produced radioactivity. In addition, by 
Varying the B,C wall thickness one can control the neutron 
energy cuu-om ease. Hetol owcetnem)/Velaw. ac aecesul tethers 
has been recently been a number of publications regarding 


ENAA utilizing boron as a thermal neutron filter (e.g. Bem 
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(source De Soete et al., 1972) 


neutron energy (in eV). 
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Figure 2.6 Cd and B cross-section variation with neutron 
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and Ryan (1980), and Stuart and Ryan (1981)). None of these 
workers apply ENAA to the analysis of geological samples 
however. In both these works, shields were fabricated by 
mixing B,C with paraffin «wax for both the outer and inner 


Sites of a SLOWPOKE reactor respectively. 


2.7.3 Experimentation and results 

Following the work of Bem and Ryan(1981) B,C shields 
were made for the outer site of the University of Alberta 
SLOWPOKE reactor. The 4 mm thick walls of the B,C saturated 
paraffin wax shields have a cut-off energy of about 10 eV 
(Rossitto et al. 1972, Stuart and Ryan 1981). To access the 
effectivness of the shields a 277.4 mg aliquot of the USGS 
standard AGV-1 was irradiated behind a shield for 10 minutes 
at maximum ttux (5.0 x 40% mn em? <s-*) an the outer sire: 
Following a decay period of 20 minutes the sample was 
counted for 600 seconds livetime at 1 cm geometry on the Win 
‘5 Ge(Li) detector. Figure 2.7 shows a segment of the 
obtained spectrum with the major peaks identified. For 
comparison, a plot of a spectrum obtained from the 
unshielded, wrradrationsof an 88.0 mg aliquot of AGV-1 1s 
shown in figure 2.8. This sample was irradiated for 4 
Minutes ate 2.5 xk 10% n cm" s 2 in Site 1 of the reactor, 
Following trradiation, except for a slightly longer decay 
time for the INAA work (30 minutes vs. 20 minutes) the 
geometry, detector, and counting periods utilized were 


identical to those employed in the ENAA of the larger 


$ 


7 


’ 


a tandyosceg 
et tabaGtades 468 SOD 
| sqnak him eWdi 5 nine. Aiod *e at nittesde itl 


: ats? mata oar 


preven Trt vovbens THA het 


et iou’ bis yol'esneatrogel® €, 
shbedde ot, 0 5 / ht amg Lee yet 16-dsGe ee palette 


sisedin 7% err ee 


lw 


G6vG4 28a '5,5) a a. 
Ve. Of giuois to’ *2o18ae 734" 
ef? -@4654 es) Pak ia 
ene “@nr $= toupblauge ?. TT 
Oma ei Dia<its "OS BNI 
oa re iii 1 ae | 
ees e | gaa 449 Pee 
Aya hi G -a@ 9 
vil ’ ! ae fie 
‘s . /Sa60! papi! &t est 


ef PPUSA 39 .QGyr'ié*on ‘ht 


| See Cae! Buy 


gi Mes a a8 f acu? rt 


Yaa 


ant); sgetain: SF 


NMwing visawksis ae 


— aSeet 


eo) 


28 asin iy Ue Rave 7 seed 
GCiow testi 103, abi rsy Cobeeiny Ehe: 
sa jol out Te sab ade oh 2 i 


43- an 6 e287 . goes 


Parr)? lal sh 803 wae worsen aa 


sti, i9Gim ety YR) oben « 


6 S¥er sbieidd ae @m3 
sage 2% la 39 Gavineel, 
Aieiae eft te spaneesoed 
to29%) 2a¢ '-e SreboageD) 
a7 
i 0.°) vf que) pee 3h 
Cc ip = Fe; eeood és veleallé 


ui inidenze® Ii2beh Pe 
- 


no echt elt ae epee Benbagdal 2 


oe sh) fete ~hand sane 
so Ve viiweibesed Sep lebdang, 
nd pilD, 1. 8-Saupl? ab-neal 
-/ me Afi eae oe A 
ye see pateaEiasa Dyas 


ene 


00 T 


1500 


1000 


DN 2 
on fe)) =a 
< E 
Le 
za N i 
o oN) S 
= = e 

(qn) E — 
oO iS = * 
= my 2 ps 
SO See 
—= be — — 
| = 
5 is 
= 
< 


750 


625 


500 


375 


250 


125 


BELO 9 + z 6elLeosg + € z 6GL9S b+ ¢& z 6eL99 + € z 6GLeS9 6 € z 


30 | v0 T sOT 20 J 10] 
S}JUNOD [eo] 


Figure 2.7 The y-ray spectrum of USGS andesite AGV-1 
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Figure 2.8 The y-ray spectrum of USGS andesite AGV-1 
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Table 2.19 Comparison of instrumental and epithermal NAA of 


AGV- 1 


INAA. of AGV-1 ENAA. of AGV-1 


Isotope keV Area Background %Error Area Background %Error 


1414 
2503 
8454 
2832 
3594 
SHAS 
US} AL 
2368 
163 


6316 


164382 


2989 


41947 


15654 


19331 


14791 


81092 
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The readers attention is drawn to the increased size of 
the **Al photopeak (produced via the reaction *°Si(n,p)?°Al1) 
following ENAA. The analysis of Si in geological samples via 
this method has, suprisingly enough (considering the 
importance of this element to the geologist) recieved little 
attention. This may be due to the short-lived nature of 
7°Al, or because $10, is routinely determined by XRF. This 
author knows of only one example in the literature where Si 
has been determined in a 'geological-like' matrix (pot 
shards) using the isotope 7°Al (Hancock 1982). 

There are however disadvantages associated with the use 
of the outer site for ENAA. Being further from the core than 
the inner sites the epithermal component of the flux is 
reduced, the gth/gepi being 58.1+1.5. A shield the size 
described contains some 8g of B.C which, due to the large 
reactivity of 7°B, probably 'draws* «the meutron flux, thus 
rendering the inner sites of little use for conventional NAA 
while carrying out ENAA in the outer site (unless flux 
monitors are simultaneously irradiated).'? At maximum, only 
two or three geological samples, each 100 mg or so in mass, 
may be simultaneously irradiated within the shielded central 
cavity (unless pelletized) thus making analysis of 
long-lived isotopes relatively expensive. 

Stuart and Ryan (1981) used an inner SLOWPOKE site and 


shield thus utilizing the greater epithermal flux component 


'2T[]) date measurement of inner site flux variations during 
the irradiation of a boron shield have not been made to 
verify this hypothesis. 
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(gth/gepi = 16.9+0.8). Their shields would also have rather 
a high reactivity as they contain some 5g B,C. Their main 
drawback is however, that samples must be placed in a 
central cavity some 6 mm diameter and 45 mm long which 
requires some delicate sample preparation, permits only one 
or two samples to be irradiated simultaneously and when 
counting, presents geometry problems. 

In an attempt to determine both uranium and thorium by 
delayed neutron counting Rosenberg (1981) irradiated samples 
mixed with and without B,C powder. From this paper came the 
idea of mixing B,C with samples for ENAA. 

It was important to determine the effectiveness of the 
ByC in reducing the thermal neutron component and to find 
empirically the optimum B,C to sample ratio. To this end 25 
mg increments of B,C from 0 to 150 mg were throughly mixed 
with 150 mg aliquots of a well homogenized coal. The samples 
were were then analyzed in accordance with the short INAA 
iradiation scheme. Figure 2.9 shows the effect on the 
activity of a number of isotopes with the increasing B,C 
content. The data are plotted aS a percentage of the 
activity of the respective isotopes in the untreated coal 
against the mass of B,C added. The error bars shown take 
into account the uncertainties in the determination of the 
photopeak areas in both the untreated and treated coals. 
Clearly it can be seen that with increasing boron content 
the activity of the major Compton producing radionuclides, 


which follow the 1/V law are supressed far more than 7°°U. 
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Figure 2.9 Eftect of BAC on the specitic activity of 


selected isotopes 
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There=seems little to be gained by imeréeasing therB,€ to 
sample ratio above 1:3 (by weight). Further experimentation 
envolved running 100 mg of shale, 100 mg shale + 30 mg B,C 
and 150 mg shale + 30 mg B,C. When analyzed it was found 
that the reduction in **Na and **Mn activity was greater 
than in the coal samples at a similar boron to sample ratio. 
Although this may be attributed to incomplete mixing in the 
coals, statistical error or that the effectivness of the B,C 
is related to its density or packing within the sample 
rather than weight, it may in part be due to the large 
percentage of organic carbon making up the coals which 
might, to a certain extent, thermalize the epithermal 
neutrons. A Similar effect due to the moderating effect of 
hydrogen in water was found by Reynolds and Mullins (1963). 
With increasing B,C (and hence C) content this might also be 
a problem. 

A die was made to pelletize mixed sample and B,C, 
producing discs 10 mm in diameter and approximately 1 mm in 
thickness for irradiation. In this way up to 20-24 samples 
(100-200 mg each) may be simultaneously irradiated in one of 
the inner SLOWPOKE sites with only about 600-700 mg B,C 
(which is much less than used by Stuart and Ryan ibid). In 
addition, flux monitors of dried Au and Mn solutions, on 
sealed filter papers, may be placed between samples thus 
allowing the tlux magnitude and 'distribution to be well 
monitored. Unfortunately, time did not permit further 


evaluation of the method. 
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The extremely well defined geometry of the pellets 
would reduce errors resulting from non-reproducible counting 
conditions. This effect is particularly important when using 
the LEPS detector where small variations in counting 
geometry are the difference between accurate and 
reproducible results and unreliable data. For a fuller 
appreciation of the importance of this effect the reader is 
refered to the paper by Hoede and Das (1981). 

Ungforgunatel yj s.sdue “to ‘the jpaucity of “tame, ithe 
preliminary work could not be followed up and consequently a 
number of questions still remain unanswered, these include:- 

bh) WS Te mixing by weight or by volume that 
controls the effectivness of the B,C ? 

2) once pelletized, will the increased packing of 
the B,C result in a more effective thermal neutron 
absorber? 

3) what are the advantage factors for the longer 
lived radioisotopes (particularly the rare earth 


elements)? 
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3. Fission Track Analysis 


3.1 Introduction 

Quantitative uranium measurements by fission track 
analysis (FTA) have been reported and the method discussed 
previously by a number of workers including Price and Walker 
(1963), Fleischer and Price (1964) and Fisher (1970). The 
method has a number of attractive features including its 
sensitivity, low cost and the variety of samples that may be 
analyzed. Full appreciation of this latter point can be 


found in Fisher (1975). 


3.2 Theoretical principles 

When a sample is irradiated in a nuclear reactor, the 
mixed neutron flux may induce (n,f) reactions (fission) in 
any heavy (Z290) nuclei present. There are however only 
three major naturally occurring fissionable isotopes, 77*U, 
Pinang ~227h SOruthe three; tonly =u" undergoest#rissron in 
a dominantly thermal flux. This is due to the large thermal 
fission cross-section of ?°25U (577b) compared to the thermal 
hisSston cross sections of “°*i and) “*-Th (lessethan Usoeand 
0.2 mb respectively, Lederer et al. (1978)). 

When ?735U fissions, two highly energetic nuclei are 
produced with masses bimodally distributed about A=95 and 
#29) (see figures? 2 75). In addition to these=two farge 
fragments, on average, two to three 'prompt' neutrons are 


also emitted. The fission fragments produced recoil in 
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essentially opposite directions to conserve both energy and 
momentum. The fragments may be recorded by placing a 
Strtable detector inPintimate contact fduesto: the short 
Rangesor, Che sEissions products) ewithethessample prior iro 
irradiation. The fission fragments create damage trails or 
EFacks) Of satomic proportions in’ the detector. After removal 
from the reactor and upon etching with a suitable reagent, 
the tracks or dislocations are enhanced to the point where 


they may be examined with an ordinary optical microscope. 


3.3 Sample Preparation and Analysis. 

Fission track studies of polished thin sections 
(Kleeman and Lovering, 1967; Hamilton, 1966 and Duke et al. 
1980) have shown that uranium is invariably heterogeneously 
distributed throughout many rock forming phases. While this 
is useful for mapping the Spatial distribution of uranium, 
samples must first be homogenized for whole rock 
quantitative analysis. Following the work of Fisher (1970) 
and Nishimura (1970) the samples of this study were first 
crushed to pass through a -120 mesh sieve and then throughly 
mixed. 

Sample aliquots 10-20 mg each, with rims and backing of 
cellulose powder, were pressed at approximately 1000 kg cm’ 
(c.14000 lbs in-?) into discs 25 mm in diameter and about 1 
mm thick. Each sample had an exposed surface area of 
approximately 50 mm?. With care, to avoid contamination, it 


is possible to pelletize three or four samples on the face 
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Lexan plastic discs were placed in contact with the 
samples and held firmly in place with scotch tape. NBL-, 
NBS- and USGS- standard materials were similarly prepared to 
the unknowns and the discs stacked in polyethelene capsules 
54 mm long with an inside diameter of 29 mm. 

All irradiations were carried out at the University of 
Alberta uSing the SLOWPOKE II reactor. The large outer site 
of the reactor, with a maximum thermal neutron flux of 5 x 
LOY Mean cis “Si eiWes mised shor the wmayority sot Teradiations,, 

After irradiation the samples were allowed to 'cool' 
por three, or, Eour days to allow the “**Na activity to decay 
to acceptable handling levels. Upon separation the Lexan 
discs were washed in water and then etched in a 6.25 N NaOH 
Soluewonat 50 Gibor.20) munutess (hasher, 197 0) 2.the solution 
was held at a constant temperature by a thermostatically 
controlled dry bath. This was found to be important as over 
etching results in damage to the polycarbonate plastic which 
makes track recognition difficult. Once etched the detectors 
were washed in distilled water and let dry. 

The track densities of the samples were calculated 
usinqean Ordinary Optical microscope, il) transmic ted sii gitar 
500x magnification (40x objective, 12.5x ocular) together 
with a graticule and point counter. 
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Where: 


geometry factor (e.1 for fagsion. tracks) 
molecular weight of uranium (238.98) 
track density observed 

1SOLODi ce abundance of “=U 072%) 
ehermad a! vux 

irradiation time 

Pission Cross Section Of - >" U(577b) 
Avogadros Number 

range of fission fragments 

density of material 

efficiency of detector in registering tracks 
weight fraction of uranium 
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(Price and Walker, 1963; Fleischer and Price, 1964). 
Determining all of the components in the above equation 
would make absolute uranium analysis by FTA extremely 
tedious. This is circumvented by applying a comparator 
technique similar in many ways to that used in INAA. This is 
achieved by ratioing the track densities of unknowns with 
those of standards irradiated and etched at the same time as 
thessamples (Fisher, 1970). This allows the cancelling of 


many of the components allowing uranium to be determined 


bye 
d(u) meen) 
C(u) = —m+:—— x C(std) ppmuU 
d(std) I(std) 
Where: 
C(u) = concentration of uranium in unknown 
C(std)= concentration of uranium in standard 
dad (u) = track density of unknown 
a(std)= track density of standard 
I (u) = Isotopic abundance of ?°°U in unknown 
I(std)= isotopic abundance of ?°°U in standard 


If we assume a constant isotopic ratio for 7°*U/?*5U, the 
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uranium content of the unknown is given by: 
d(u) 
C(u) = ————— x C(std) 
d(std) 

The track densities were calculated by counting ten 
randomly selected, equal areas and taking their mean value. 
The standard error of the means of the known and unknown 
Samples were calculated and the total error or uncertainty 
in the uranium being the square of the sum of squares of 
these two errors. 

The total uncertainty is largely controlled by the 
known and unknown standard deviations which in turn reflect 
the spread of the data points about their respective means. 
Consequently, assuming a statistically significant number of 
tracks produced (which is proportional to the integrated 
neutron flux) and counted, the precision of the method is 
largely controlled by how well homogenized the samples are. 
Where results with a precision consistently <+5 % are 
required, standard glasses are available for comparison 
(supplied by the U.S. National Bureau Of Standards in four 


concentrations : 461 ppm, 37.4 ppm, 0.82 ppm and 0.072 ppm 


Urs 2 


13 There is however some doubt regarding the homogeneity of 
the uranium distribution in the 0.072 ppm U standard 
(Fisher, 1975 p.295) 
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3.4 Results 

Twenty Exshaw Shale samples were prepared and analyzed 
by FTA. As can be seen from a comparison of the results of 
the various methods of uranium analysis given in table 7.1 
PTA 1S sufficiently accurate for’ the routine analysis of 
uranium in geological materials. The method does not suffer 
any major matrix effects, and is very sensitive. The results 
are however, rather imprecise. This is largely a consequence 
of heterogeneities found in the NBL-81A (20 ppm U) 
comparator standard used, as exemplified by fission track 
"stars' in the detectors. Similar stars were observed in 
some of the Exshaw Shales, particularly those with higher 
than average uranium concentrations. Use of the NBS-glasses 
would avoid this problem. 

Theoretically, if the material used to register the 
fission tracks has an extremely low U concentration then DNC 
is one of the most sensitive metods available for uranium 
analysis. The method requires very little sample and may be 
applied over 5 to 6 orders of magnitude of uranium 


concentration. 
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4, Delayed Neutron Counting. 


a, /)introduction 

Approximately 110 samples (including replicates) were 
analyzed for uranium by delayed neutron counting (DNC) at 
McMaster Nuclear Reactor (MNR) and later at the SLOWPOKE 
reactor. Only a short review of the principles will be given 
here and the reader is referred to Amiel (1962), Gale (1967) 
and Boulanger et al. (1975) for a more complete discussion 


of the theory and practice of DNC. 


4,2 Theory 

As noted when discussing the principles behind fission 
track lanalysis-=utiwas pointed ‘out that faissionvote--Usis 
selectively induced by thermal neutrons when a sample is 
irradiated in a nuclear reactor. The fission fragments 
produced are unstable and decay with relatively short 
half-lives (maximum of approximately 55 s) with the emission 
of B--particles. In some of these decays (1.58% Keepin et 
aln9(1957)) this is followed by the immediate emissionsotve 
'delayed' neutron. These delayed neutrons, unlike the prompt 
neutrons emitted during the fission event, continue to be 
emitted after irradiation. The number of delayed neutrons 
produced 1s proportional to the amount of? “U presen: and 
when counted can provide a rapid, low cost, precise, 
accurate and extremely sensitive means of uranium analysis. 


There are however some light radionuclides, produced by fast 
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neutron reactions, which also decay with the emission of 
delayed neutrons. Of these °Li (T1/2=0.17 s) and ‘7N 
(T1/2=4,.17 s) produced in the reactions *Be(n,p)°?Li and 
“"O\n,p)"7N are the most important. The possible 
interference due to these two isotopes is discussed later. 

ASeit is actually ~2Flethateis determinedevias DNC, for 
quantatitive total uranium analyses, it must be assumed that 
Bhes4* *U/2*"U-isotopiceratiosis «constants forneunate!y,owhen 
analyzing geological samples, this is nearly always the 
Gasca & 

Amie leibid ecaleulatedsthat theresismlitcie tombe 
gained by irradiating a sample for more than one minute. The 
reasons for this are twofold. Firstly, due to the short 
half-lives of the neutron emitting fission products, the 
Signal—to-background natio 1S snot increased «significantly by 
a longer irradiation. Secondly, the induced (n,y) activity 
(largely due to 7*Al, T1/2 = 2.24 min) may interfere when 
counting and also presents a potential radiation hazard. 

A counting period similar in length to the irradiation 
time is recommended as it allows the majority of the delayed 
neutron emissions to occur and be detected. Counting for 
shorter periods of time (as compared to the irradiation 
time) increases the sample throughput but reduces the 
Sensitivity of the method. Too long a counting time on the 


other hand reduces the precision by increasing the 


‘4 There is only one well documented variation from this 
Situation, at Oklo in Gabon where *°*°U is depleted due to 
its natural fission (see Bouzigues et al., 1975). 
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registering slow neutrons are enriched BF; proportional gas 
counters. The meutrons interact with’ the enriched %°B to 
producejalpha particles and "hi nucle? via the reaction 
"p(n, a) Li (thermal cross-section of 3840 barns). Being 
positively charged, alpha particles are readily detected 
when a high voltage is applied across the detector. Due to 
the extremely large thermal cross-section of the above 
reaction (and small fast neutron cross-section) maximum 
detector efficiency only occurs when the higher energy 
delayed neutrons have been slowed down to thermal energies. 
This is accomplished by placing the irradiated sample in the 
centre of an array of BF; tubes (8 at MNR) set in paraffin 
wax or polyethylene. The latter materials, being rich in 
hydrogen, slow down the delayed neutrons by elastic 
scattering, thus increasing their chance of detection. 

On account of the large reaction cross-Section and 
availibility of enriched boron, the BF; detectors are 
comparatively speaking, highly efficient and well-suited for 
this application. The voltage pulses generated during the 
detection of the neutrons are passed through a pre-amplifier 
and amplifier to a single-channel analyzer (SCA) where they 
ane accumulated. The gross ‘counts at the end of the count 
interval may be read directly or output to a line printer. 
Once the system is calibrated against a number of uranium 
standards of varying concenratations (of known isotopic 


abundance) samples may be run on a routine basis. 
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Towards the end of this study a DNC system was built at 


the University of Alberta SLOWPOKE facility for uranium 


The SLOWPOKE counting system consists o 
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errey <5 cMS in diameter and embedded in paraffin wax. A 
jead cylinder (8 cms diameter, 25 cms long) was placed 


coaxially in the middle of the BF; tube array to attenuate 
the gamma photons due predominantly to the short-lived 
radioisotope **Al (T1/2=2.24 min). The whole counting system 


was Surrounded by a further 


an) 


cms of borated paraffin wax in 
an attempt to reduce the external background component. 

The BF; tubes were connected in parallel and operated 
at a common high voltage of 2200 V supplied by a Tennelec 
preamplifer to a Tennelec TC216 linear amplifer and SCA and 
on to a TC534 counter and timer. Raw data were printed out 
on a decwriter via a Tennelec TC588 buffer printout 
control/interface. 

The discriminator bias on the counter/timer was 
adjusted with the aid of a Canberra Series 30 MCA so that 
all of the electronic noise and low amplitude pulses 
produced by the interaction of gamma-rays and the detector 


tubes were eliminated with only a small reduction in the 


neutron detector efficiency. 
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4.3 DNC results 
Preliminary measurements, using an irradiation : decay 
GoOuntesohemesore20e10e20sS.aatya  fiuvwore | xaos neem = 
Se edaverta Sensitivityrot o0,it iad ctseuga® fRaturalsuraniun, 
Decision and detection limits (Currie;1968) for uranium are 
50 and 120 ppb respectively in a 5g sample. 

Analysis of a liquid atomic absorption thorium standard 
under identical conditions gave a sensitivity for thorium of 
CeooHUsO Sect Sedass that Book on e006 pe thee of tirand um. 
As expected, this ratio is in good agreement with the 
Sensitivity ratio of 1.17%+0.05% determined by AECL 
Commercial Products (see Boulanger et al. 1975) who also 
used a SLOWPOKE reactor (but different irradiation, decay 
and count times). 

Using the Canadian Radioactive Reference Standard BL-1 
(U=220410 ppm; Th=15+1 ppm) as the calibration standard, the 
accuracy of the method was assessed by analyzing a number of 
reference standards, the results of which are given in table 
ara le 


Table 4.1 DNC uranium results using the SLOWPOKE reactor 


Sample Th/U U( ppm) Corrected Recommended 
ID Ratio (+1 sigma) t U (ppm) * value 
USGS-W1 Bie. 0 OR Grl0s. Oi OPS 9AO ROS O58 
Di 2.0 AOE Nee Bono tales 41+2 
Dizi ZnO AAO ease 2 CUS nal Pal 41+2 
NBS 1632a 2.4 10.920 .4 10.6+0.4 Toe OF 
NBSei635a 3.5 Lo Soe 2 (e282 0Rel 2 eo Oral 2 


+ results of a single analysis of each standard. 
* correcting for the known Th content of the sample. 
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As can be seen, even without correction for the thorium 
content of each sample, the uranium results are quite 
accurate. Following correction, all results, except NBS 
1632a, lie within one standard deviation of the recommended 
uranium values. 

The thorium interference need only be considered 
Significant when analyzing samples in which the Th/U ratio 
is greater than about 5. For such samples, determination of 
the thorium content, and hence the interference to uranium, 
should be made. One possible method is to re-irradiate each 
sample in a predominantly fast neutron flux (see Amiel ibid. 
and Gale ibid.). This effect may be seen in figure 4.1 where 
the uranium results obtained by delayed neutron counting at 
McMaster Nuclear Reactor are plotted against those obtained 
at the SLOWPOKE reactor. In general, the values obtained at 
SLOWPOKE are consistently higher than those from MNR due to 
the greater gth/gfast at SLOWPOKE (approximately 4.6) 
compared to the ratio at MNR (14.6 pers. comm. E. Hoffman). 

Millard and Keaton (1982) throughly investigated the 
factors affecting the error in the DNC technique for U and 
Th analyses and concluded that counting statistics are the 
MajOr Source of uncertainty. Consequently, under the present 
scheme a concentration of 0.4 ppm uranium in a 5g sample may 
be determined with a precision of better than +10% (one 
sigma). 

Replicate analyses of five aliquots of the same sample 


analyzed using the SLOWPOKE DNC system gave a mean value of 


_ 


Siu lage 4 4.02 


| | eee 
1 wabetile @02 Jol oeltae 
. Ca aa ft 

4h fats evra lis eal ivartes ans 
LytiisamedsT Ste I2 \cTanteat basins am nds 


etre /ave®..\ 0142 ofa Bae 
Hiei $6 GALI) Gir AOTw ott 


: Sur witu: 
42nd Drlt wis an 
: Ags salt b1NO? 524)\ a ekunae ast; 78 i 
py | Ate evi dea gp o-sa 4 
PAR @GItnil&e? ~ @ ric Hots 
okt aw | , : ( canes 24° 
Ae LS | 7) 
ia. wisaei ( 5 Pimse n 
mur ese SS) - ee te re La 
tw henidé wife cel /b4S8 09%) 9 
=" as 4 ai 


hn iat } 


en! frr.!s 

Rvit 

ees +h 
See te oN: 


CHG aby @ bem 


eiqans oan. 04316 Waipi ia as® \anapeyle ‘2 
fo une dnee @ aver findeve Um AACA of , a - 


fq@¢ hers yi tel “aud 


a 


> wld 489902 io hol cies GONe/ Geaderegey ad 


a Le 


V4 ea ei te 


On a IG om 


PAT Ri4Ago- 
OD Ihevpaena 


misnsw # a! 


' 


hh 


a 


pone ees pri sugtss e70)h03 


fi | - 


a: Me os 


aes elgese 


oo. eu h iprey bites ts cis. SEZ PIAS wprhende elt 
LMienslawe once ihe ante yraobté: 


rie Solid 


A 


= 


beaker 4 


= sutdu anise 
17? Aa wid 10d 
he ad btu : 


> suan{esewig ev) wig i 


P7987 wae ln ae 
cisgs7 DHOAMDIE 443 aa’ 
i@rbscTrhs> @36 seoaware | 
AVE neq Tenet fa 7 

of dh OP. HEF OT bai egm it 


anges) San. peal liz 


ieineLomc> Gis wee lane 


*iierne tate he apIeae wie 


2 Ml SORIIEISS @ Stee aa 


150. 


-DNC) 
o 


U in ppm (SLOWPOKE 


—s 
Oo 
© 


~ 
os 


on 
2 


No 
oI 


81 


45 


GE 


9 3 Sigma Errors 


ZS. 90) 12 10032 oe LOU: 


U in ppm (MNR - DNC) 


Figure 4.1 Comparison of SLOWPOKE and McMaster reactors DNC 


uranium results 
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five samples analyzed at MNR gave a mean value of 27.01+0.15 
Pompe OF oo, FN=Spe seartere=0 (n10%jfaClecarly then, as 4 
method of uranium analysis, DNC is extremely reproducible in 
Neacuae, 

With the present irradiation:decay:count scheme it is 
possible to analyze some 60 samples per hour which, with the 
operating constraint of 4 hours at maximum flux at the 
University of Alberta SLOWPOKE reactor allows for some 240 
samples to be analyzed per day. 

A FORTRAN IV programme was written by the author which 
given the sample identity, mass and raw counts for each 
sample and calculates the mass and concentration of uranium 
together with the total error of the analysis (based on 
background removal, counting statistics and division by the 
standard sensitivity). 

For an excellent coverage of the potential 
interferences to the analysis of uranium by DNC the reader 
is referred to Kunzendorf and Lgvborg (1981). They conclude 
that the method is fast and reliable, and well-suited for 
the analysis of all kinds of geological samples. In their 
work they found errors due to variations in timing together 
with sample preparation negligible and enhancement of 
Fission due to “=’Pul or 27°7U in normal geological samples 
Unlikely. The typical concentrations of B and Cd in 
geological samples was found to have a negligible effect on 


fhe neutron flux. Absorption of ‘thermal neutrons was only 
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found to be a problem when the rare earth elements Sm and Gd 
Were Present at concentrations levels of the order of i000 


ppm and then corrections were necessary. 
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5. Gamma-ray Spectrometry. 


5.1 Introduction 

Traditionally when analyzing rocks, sediments and soils 
for uranium and thorium by gamma-spectrometry, only the 
higher energy gamma emissions of the daughter products are 
usually utilized as they are more interference-free. 
Typically uranium has been quantified indirectly using the 
V./6 MeV “Gamma-ray of °"*“Bi > and@thorium indirectly via the 
2.62 MeV gamma emission of ?°*Tl (see figure 5.1 for the 
7, = Th and 2-1) decay chains). The wteasonstor tars sare 
essentrallyithree-toOlde Firstweot alle 29") andere ih 
themselves decay by a-emission with insignificant gamma 
production.'*® Secondly, few of their immediate daughters 
which do decay by gamma emission emit sufficiently intense 
gamma radiation for quantification. Finally, during 
exploration for hidden mineralization only the higher energy 
gamma emissions may be used as the lower energy photons are 
totally absorbed by a thin soil cover. Even if this 
absorption did not occur, the low energy spectrum would ride 
on a significant low energy background which would decrease 
the sensitivity of the method appreciably. 

Consequently, when determining U and Th via their 
daughter products the assumption is usually made that the 


15 2351 does decay directly emitting a low energy gamma 
photon at 0.18 MeV which has been used to quantify U in the 
laboratory (Bunker and Bush, 1966). The method is limited to 
Samoles with a Th/U ratio <10. Above this) the interterence 
due to ??2%Ac (a ?2?Th daughter product) is too great for 
accurate U determination. 
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parent radionuclides are in secular equilibrium with their 
daughters. In secular equlibrium the number of atoms of each 
daughter being produced is equal to the number of atoms of 
the same daughter decaying. This requires that the system 
remain closed to gains and/or losses of radioisotopes of the 
decay chain under consideration. When such conditions are 
not met disequilibrium results. Some eight or nine 
half-lives of the 'offending' radioisotope are then 
necessary for equilibrium to be regained. 

Due to the relatively short half-lives of the decay 
Products einsthe ==47h decay chain. disequilibrium in this 
chain is rarely a problem.'* However, many of the daughter 
Broducts inthe "UU decay, chain are sutftiecrently long-lived 
and chemically distinct that selective mobilization, by 
leaching for example, can result in loss of equilibrium that 
is not rapidly regained (see Levinson and Coetyee, 1978., 
Szogny and Kish, 1978, and sRosholt,. 1959). 

Only disequilibrium due to the diffusion of *??Rn (part 
of the 72*U decay chain) is of immediate concern in the 
laboratory (see under "Sample Preparation"). A comparison of 
the uranium results obtained by gamma-spectrometry with 
those from INAA or DNC (where uranium is determined 
directly) will identify samples in which disequilibrium is a 
problem. 


i¢222'Pashas the Longest half-life of sche daughter products 
in the 22?Th decay scheme (11/2 = 6.7 yr). Hence, >99% 
equilibrium is attained after 46 years (approximately seven 
half-lives). 
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Potassium may be determined via its monoenergetic gamma 
emission at 1.46 MeV due to the decay of ‘°K following 
electron capture. AS potassium is determined directly by 


gamma-spectrometry disequilibrium is not a problem. 


5.2 Sample Preparation 

Samples to be analyzed were washed briefly with tap 
water to remove any traces of drilling mud. Then, after 
drying, they were crushed in a tungsten-carbide swing mill 
to approximately -120 mesh. They were then vacuum dried for 
24 hours and weighed into tared plastic vials (5.4 cm long 
with an inside diameter of 1.4 cm) to a height of 4 cm. This 
was found to be important for two reasons. Firstly, a 
non-reproducible counting geometry will result in 
non-reproducible results. Secondly, and having the same 
effect, the efficiency of the NaI(Tl1) well detector utilized 
was found to decrease rapidly towards the top of the well as 
the near 4% geometry was lost. To avoid later settling the 
vials were tapped to the above height during packing. The 
mass of crushed rock was measured to the nearest milligram 
and recorded. 

The plastic vials were then hermetically sealed and 
left for some thirty days before being counted. This is 
necessary as it allows the ***U decay chain to regain 
Secular equilibrium as; when crushed, “the gaseous © °" Rn 
daughter is lost. Consequently, °*“Bi which occurs even 


later in the decay chain (see figure 5.1) may only be used 
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to quantify the uranium once equilibrium has been regained. 
In this work a’minimum of seven half-lives of *22Rn (T1/2 = 
3.825 days) were allowed for the gas and its daughter 
products to build up and regain >99% equilibrium. However, 
Ritchie and McHenry (1973) found from their work on uranium 
analysis by gamma spectrometry that: 

"measured gamma-ray equilibrium was achieved 

im 40-15 days.” (Ritchie and McHenry p57 7 

Log 208 
Radcte 2209 (Ji 26=) S5es Secs) ocecurs in theo?" Th decay 


chain, however, due to its short half-life, equilibrium is 


regained approximately eight minutes after sealing. 


5.3 The Counting System. 

Samples were counted using a 3" x 3" well-type Bicron 
Nal(Tl) crystal detector with a measured resolution of 7.65% 
at the '*’7Cs 0.662 MeV photopeak. The crystal was mounted on 
a photomultiplier tube which was in turn connected to an 
Ortec 276 photomultiplier base and preamplifier. A lead 
castle some 10 cm in thickness housed the detector system. 
The 1 kV detector bias was supplied by an Ortec 456 high 
voltage power supply. The preamp signals were fed to a 
Canberra Series 80 multichannel analyzer (MCA) where they 
were amplified and assigned to one of 1024 channels. Spectra 
thus obtained were stored on floppy-discs using a Canberra 
Model 8660 Series 80 Floppy-Disc System for later numerical 


analysis. 
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Any drift in the system was monitored using an Ortec 


419 precision pulse generator. 


5.4 Analysis 

As noted earlier U and Th have generally been 
determined via the high energy gamma emissions of 7'‘Bi and 
Tl (Adams et al... 1956). Although conventent for field 
analysis this method has two main drawbacks for routine 
laboratory work. These are: 

a) large samples and 

b) long counting times 
are necessary when measuring low levels of U and Th (0-50 
ppm). Figures 5.2 and 5.3 show typical spectra of uranium 
and thorium standards over the range of 0 to 3 MeV. From 
them it can be seen that the low energy regions of the 
spectra (0 to 0.5 MeV) contain the greatest amount of 
potential information. This is duewin part to the greater 
intensity of the low energy peaks, together with the 
increased efficiency of the detector at about 0.1 MeV. Both 
factors result in an increased countrate and better counting 
Statistics. 

After a series of preliminary studies at the 
concentrations of uranium expected in the samples of this 
study (0-100 ppm) it was decided that ‘conventional’ gamma 
spectrometry would not be sufficiently sensitive or rapid 
for the routine analysis of a large number of samples. It 


was therefore decided to investigate the lower part of the 
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Figure 5.2 Uranium spectrum (0 to 3000 kev). 
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Figure 5.3 Thorium spectrum (0 to 3000 kev). 
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spectrum. 

In an attempt to utilize the mass of information 
contained in the lower energy part of the spectrum (see 
Hurley, 1954, 1956; Bunker and Bush, 1966: Culbert and 
Leighton, 1981) the method of Culbert and Leighton was 
adopted for U and Th meaSurement. Following their work the 
amplifier gain was adjusted so that only the lower energy 
part of the gamma spectrum (0-400 keV)'’ was assigned to the 
1024 channels of the MCA. Two windows, or regions of 
interest (53-112 keV and 208-269 keV) were set up to 
encompass the 80 keV photopeak of *'*Pb and 235 kev 
photopeak=ot “Pb. (daughters iat 72") and) @e 2th 
respectively). Figure 5.4 shows spectra of the uranium and 
thorium standards superimposed, with the photopeaks and 
windows identified. As can be seen, neither window is 
interference-free and consequently stripping coefficients or 
ratios must be determined. Knowing the mass and U and Th 
concentrations of two standards, the stripping coefficients 
for U and Th were calculated after repeatedly counting the 
background and standards. For the equipment and operating 


conditions employed, the coefficients were calculated to 


bes— 
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GPU s= 917 396>* 010 ectcrs = d- ppm =U 
Cit = W604 x 20-“~cts sg pom ath 
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17 While discussing low energy gamma spectrometry y-ray 
energies will be given in kilo electron volts (keV) for 
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Figure 5.4 The low energy parts of the uranium and thorium 
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Onee thes background counts havesbeem subtracted, the 
uranium and thorium concentrations of the sample may be 
calculated from the following equations:- 


(N}2S2Te— Ne. Str) 


U (ppm) = 
i Sot ae oT RBS OU aS LT) 
NeeS2u a Ne eSiu) 
Th (ppm) = 
meCS2teS 1 Ses iUes on) 
where:- 
N;, = net counts per second in window 1 
Nz = net counts per second in window 2 
m = mass of sample in grams 
SnX = Stripping coefficient for element X in window n 


Following this method, the detection limits (Currie, 
1968) in a thorium- and uranium-free matrix are 0.63 ppm U 
and 1.57 ppm Th respectively, for an 8-9 gram sample counted 
for 10,000 seconds. 

Such short counting times, together with the above 
sensitivity, are a great improvement on conventional 
gamma-spectrometry. In order to increase the sensitivity 
further, either the sample size must be increased and/or the 
background decreased. Much of the background in window 1 is 
due to 75 keV fluorescent Pb Ka X-rays produced by the 
interaction of cosmic radiation and high energy background 
gamma-rays with the lead shield. In an attempt to reduce 
this component, a copper shield (1.25 mm thick) was made to 
Surround the detector to absorb the Pb Ka Z-rays. A 
background reduction of approximately 33% was found in 


window 1 when the shield was used. With a thicker copper 
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shield this reduction could be further increased. 

One major problem associated with this method is the 
contribution to windows 1 and 2 by bremsstralung X-rays. 
These X-rays result from the kinetic energy loss during 
interaction of energetic £°-particles (produced during the 
decay of *°K) with the sample matrix. To examine the extent 
of this interference, a 6g sample of reagent grade KCl was 
prepared and run in accordance with the above method. From 
EbUSemt was ‘calculated that wz Kyl. 205% Ke0) an ened 
sample results in an apparent uranium and thorium 
concentration of 0.75 ppm and 0.34 ppm respectively. 
Typically shales contain between 0 and 4% K,20, therefore 
uncorrected U and Th determinations could be high by as much 
as 3 ppm U and 1.5 ppm Th. This does not represent a 
Significant error when measuring U and Th levels >50 ppm. 
However, it should not be ignored when meaSuring lower 
levels, particularly in the 1-10 ppm range. 

As the samples of the study could lie in this range, it 
was clear that a method to correct for the potassium content 
needed to be developed. Culbert and Leighton (ibid.) suggest 
that the induced fluorescent photopeak of iodine (29.8 kev) 
in the NaI crystal could be used to monitor the potassium 
interference. However, this procedure was found inadequate 
by the present author. For the purpose of this study, it was 
decided that the low energy gamma-spectrometry method be 
modified to include the analysis of K and hence correct for 


the interference. Alternatively, the potassium content of 
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the sample may be determined by an independent method such 
as &-ray fluorescence spectrometry and a correction made in 
conjunction with the above empirical measurements. 

In Order tO measure the potassium content, the gain of 
the amplifer was altered so that the range 0 to 1600 keV was 
assigned to the 1024 channels of the MCA. A third window 
(See = oures®5.5, 5976 and 5.7) was Set up to encompass ene 
14600 KeVephotopeak of “°K im addition to the Ueand Th 
windows discussed above. Extended counts (100,000 seconds) 
with an empty vial placed in the detector well were 
repeatedly taken so as to reduce the error in the 
determination of the background component of the windows. As 
in the case of the low energy gamma-spectrometry, stripping 
coefficients and equations for calculating the U, Th and K 
contents of samples had to be determined. Repeated 
measurements on a series of U, Th and K standards were 
preformed. Analysis of these measurements yielded the 


following Stripping coetficients:— 


Slee lg x Ore lSe Ss eG apple .U 
C20 = ie Ooo x LO Ch seS so) Opi ae a 
S3Ue= Se oa LOS Se rSeSa Gg. poms au 
Git =o 12054 x L0e4etse se’ go ppm Arh 
C20 e= Ooo x LOn Cc Uses gO) | Dis wld 
Gor] 6.447 x 10Setsisy Gg “pon ein 
Cike=re2 Sa xe Un ecm Se OU wie 
G2ke= e220 KURO Ct ores. “Gal, sh 
Sokt=03, 020.58 (0 SChapoe "GO, Vike 


The U, Th and K content of samples may be determined from 


the equations on page 97. 


Where three 'pure' standards are not available the 


stripping coefficients may be calculated from three mixed 
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Figure 5.5 Uranium spectrum (0 to 1500 kev) 
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Figure 5.6 Thorium spectrum (0 to 1500 kev) 
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Figure 5.7 Potassium spectrum (0 to 1500 kev) 
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samples/standards of well known U, Th and K content. 

To simplify the numerical analysis of the results, an 
interactive FORTRAN IV programme was written by the the 
author. Given the counting time, integrated counts in 
windows 1, 2 and 3, plus the sample mass, the programme 
Calculates the U, Th and h concentrations, errors, ratios: 


etc. for each sample. 


5.5 Results 

To examine the accuracy and precision of the low energy 
gamma-spectroscopy method for uranium and thorium analyses 
the reference standards BL-1, DL-1 and NBS 1633a (coal 
fly-ash) were prepared and analyzed a number of times over a 
three month period. The recommended and meaSured uranium and 
thorium values are given in table 5.1. 

Uranium determinations for the radioactive reference 
standard BL-1 are in each case well within the recommended 
limits. For the same sample however, the accuracy of the 
thorium determinations varies and in all cases the one sigma 
uncertainty 250% (usually much greater)? The targe 
uncertainty is due to the unfavourable Th/U ratio (0.07) of 
the standard. The method is most precise for both thorium 
and uranium determinations when the Th/U ratio lies between 
0.5 and 2.5 (this being well exemplified by the analyses of 
Diedei. 

When counted for 10,000 seconds, uranium and thorium 


concentrations down to 1 and 3 ppm respectively may be 
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Table 5.1 Uranium and thorium results of low energy 


a) 
b) 


Cc) 


rec. 


a) 


b) 


a) 


b) 


gamma-spectrometry 


Canadian Radioactive Reference Ore BL-l 


Uranium’ (ppm) Thorium (ppm) T count (s) 
PACE) Panis) Oeste Geo 10 ,000 

SHEA hse Fi Ded. Sew), 1000 

224228 Wey ae FLO bs Al, 600 
2202 oO Stel: 


Canadian Radioactive Reference Ore Dirl 


Uranium (ppm) Thorium (ppm) T count (s) 
Boag etel.S eeaget aa T1\0) 10,000 
20.56 2 1.9 BOG ut 2.2 7200 

Agi 2 (Skee ae, 28) 


NBS Flyash 1633a. 


Uranium (ppm) Thorium (ppm) TeCounize (S) 
132952 0.6 ZaSOut sO. 36,000 
124 220 7 2420t. 0 10,000 
LO ete Oia: 245] 0.8 


(all values plus or minus one standard deviation) 
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determined with a precision better than +10% in a thorium- 
and uranium-free matrix respectively. 

The variable counting times necessary for analyses show 
that accurate and precise results may be obtained with short 
(e.g. 10 minute) counting intervals when the uranium and/or 
thorium concentrations are sufficiently high (>100ppm). 
Consequently, as noted by Culbert and Leighton ibid. the 
procedure may be used to rapidly 'screen' samples in an 
exploration programme. 

The thorium results obtained for NBS 1633a are in good 
agreement with the recommended values, while the uranium 
values are in each case high. This is thought to be due to 
the potassium content of the sample. It was calculated that 
with the recommended 1.88% K of NBS 1633a with the sample 
mass used (7.822g) would produce the equivalent uranium 
eoncentrati1onemof Fis.5°+ 0.tppm. When “correctedutormeth iis 
interference, one of the two uranium analyses falls within 
the recommended limits. The reason the other result is some 
10% high is unclear, but it may be due to temporal 
background variations. Both the uranium and thorium results 
obtained for the Canadian Radioactive Reference Ore DL-1 are 
in excellent agreement with the recommended values, before 
and after removal of potassium equivalent uranium. 

Comparison of the uranium results determined via the 
LEGS and DNC (MNR) methods (see chapter VII) show the 
agreement to be excellent. This agreement is also shown 


Graphically in figure 7.5. (Erom this figuoe wit can (a lsovbe 
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seen that in many instances, at uranium levels below 40ppm, 
the results obtained by gamma-spectrometry are slightly 
greater than those values obtained from the DNC method. This 
is largely due to the potassium interference to the uranium 
window. To assess the accuracy and precision of the low 
energy gamma-spectrometry with potassium being quantified in 
addition to uranium and thorium, NBS 1633a was analysed, the 
results of which are’ given in table 5.2. Although the 
uranium results are in each case within the limits of the 
recommended value, they are consistently high by a factor 
1.053 + 3.41% (N=4), scatter = 2.33%. Using the inverse of 
this value as a correction factor and applying it to the 
uranium analyses, a 'corrected' weighted mean of 10.20 + 
Oe3s9pom UW (N=4)> Scatter = 2.534% results. Thus the 
eorrection factor 0.95 2 6.241% has been applied to the 
uranium results determined via the 3 window Nal method. 
Similarly, from the same four analyses the thorium 
results were found to be consistently in error by a factor 


of 1.08 + 1.83%. The inverse value, 0.93 + 1.83 (N=4) 


Scatter 2.83% was used to “correct” “the results it us 
difficult to assess the accuracy of the thorium results as 
the ‘real' thorium contents of the Exshaw Shale samples are 
UnCenia IN. 

Analysis of the potassium results of NBS 1633a show the 
values obtained in this work to be greater than the 
recommended value by a factor 1.063. Using these results the 


correction factor 0.941 + 4.10% (N=4), scatter = 9.77% was 
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Table 5.2 Low energy gamma-Spectrometry for the analysis of 


KO 


uranium, thorium and potassium in NBS Flyash 1633a 
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calculated and applied to the NaI values. Once applied the 
agreement between the results of the NaI and XRF methods are 
Generally iwithin thescaiculated wimuts (see wtable. 7,1). 

In summary therefore, both gamma-spectrometry methods 
provide a relatively cheap (as neither a nuclear reactor, 
semiconductor detector nor sophisticated data reduction 
system is necessary) and rapid means of uranium and thorium 
analysis. The LEGS method is both accurate and precise for 
both elements at normal Th/U ratios, and may be used for 
rapid screening of samples and more precise analyses. The 
three window gamma-Spectrometry method offers the advantage 
of being able to analyze simultaneously for all three of the 
major heat producing elements at little or no extra expense. 
The precision or uncertainty in such meaSurements is largely 
controlled by the counting time and relative concentrations 
of uranium, thorium and potassium of each sample. Two 
sources of random or unsystematic error may manifest 
themselves when analyzing for uranium by either method. 
Firstly, a disequilibrium break in the 2-21) "decay  chadm sau, 
or above ??*Ra, will give spurious results. This would 
become apparent if the sample was re-analyzed by INAA or DNC 
where uranium is measured directly. Secondly, if the samples 
to be analyzed have a high average atomic number (e.g. 
mineralized or panned heavy concentrate samples) attenuation 
and/or self-absorption of the low energy radiation of window 
1 (7324Th, 2'4Pb plus Bi conversion X-rays) will occur. 


Self-absorption would be readily apparent from the high 
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count rates observed and the higher energy, more 
conventional gamma-spectrometry method may be used. Sample 
density, and hence mass, (as vials are packed to the same 
height as closely as possible) may be used to identify those 
Samples in which matrix absorption may be a problem and 


those samples analyzed by an alternative method. 
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6. X-Ray Fluorescence 


Ceieintroductron 

All X-ray fluorescence (XRF) analyses were carried out 
commercially at the University of Nottingham, England using 
a fully automated Philips PW 1400 wavelength dispersive 
SpeceLometer wand “rhodium tube, vA full deseription of “che 
methodology and corrections employed may be found in Leake 
Gora i965) fand tHarvey@#andraAtkin (1962)0 

For accurate major oxide determinations sample aliquots 
of 1-2g were fused to make glass beads. The trace element 
concentrations (and in some instances the Si and K contents) 
of the samples were determined via pressed pellets. 
Approximately 7g of < 200 um of homogenized material was 


required to produce a disc of 'infinite' thickness, needed 


for the short Rh wavelengths. 


6.2 Analyses 
Five samples (of various U content) were analyzed for 
Enesmajor Oxides to 3- 
ue classify the sedimentary pelites 
De see if there was any obvious correlation between 
U, and Th, and any of the major oxides. 
on compare the results with those of INAA. 
The same five samples, plus another 35 were analyzed for the 
trace elements Cu, Pb, Zn, Th and U (and in addition, the 35 


samples were also analyzed for Si0O> and KO). The WU, Th-and 
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K contents were determined for direct comparison with the 
other techniques employed, while Cu, Pb and Zn were 
determined because they are difficult to analyze for via 
INAA at trace levels. 

The reported approximate 2 sigma detection limits for 
the XRF analysis, based on a variety of rock types for a 
number of trace elements is given in table 6.1 (written 
communication B. Atkin). It can be seen that the detection 
limits for U and Th by XRF are not as low as those of DNC, 
INAA or FTA. However, Galson et al. (1982) have been able to 
lower the 2 sigma detection limits to 0.2 ppm U, 1.0 ppm Th 
and 0.001 % K with better precision than previously 
obtained, but still uSing a single pressed pellet. These 
results are obtained using analytical counting times of 
4800s, 480s and 72s for U, Th and K respectively. Increasing 
the counting times by a factor N would lower these detection 


limits by the square-roor of N: 
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Jable 6.) Detection 1Timits for X-Ray fluorescense 


Element Lower limit of detection 
Ag* 10) em, 
As Lepomis 
Ba 8.3 Dum. 
Bi i Ospem. 
Ce WO) em. 
Co 2.5) PEM. 
Ge 2.0) DEN. 
Cu 0.8 ppm. 
Mo 0.9 ppm. 
Nb 1.2 ppm. 
Ni he a jajein, 
Pb Ae 8/2 116 
Rb 0.7 ppm. 
S 8.0 ppm. 
Sb* i342 pom. 
Srey 2.0 ppm. 
Sr O29) pan. 
Th 0.6 ppm. 
U 1.4 ppm. 
V EAS eje.i0e 
W 1.6 pr. 
mn 3.0 ppm. 
Zc Siolh (eis ite 


(* - Cr X-ray tube; all other elements Rh X-ray 
excitation) 


(written communication B. Atkin, 1982) 
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7. Evaluation of the various methods for uranium analysis 

Which sf ethe various «methods. for uranium analysis 
employed in this study is the 'best'? Such a question does 
not haveva straightforward answer as the term “best’ is 
ambiguous. Does it mean the most precise, most sensitive, 
cheapest, quickest or simplest method, to name a few 
possible interpretations? Clearly, the 'best' method depends 
upon the requirements of the analysis and the availability 
of the various techniques. 

For example, in the case of a geochemical exploration 
programme for uranium, where there are many thousands of 
samples to be analyzed solely for uranium, then the easily 
automated delayed neutron counting method should be employed 
on account of its speed, senSitivity, lack of interferences 
and low cost. However, aS 1S apparent from the earlier 
discussion on the technique, access to a nuclear reactor is 
essential. If only very small samples (<20 mg) are available 
fou anmalysis* jand/or ‘thesuranium convent fofethe sample is 
expected to be very low (say <500 ppb) for example in 
ultramafic or tholeiitic rocks, then the extremely sensitive 
fission track method should be used (Fisher, 1970,1975). 
Once again however, a reactor is necessary to bombard the 
samples. 

For samples where more than U and possibly Th are to be 
determined, DNC and FTA will not be able to satisfy all the 


requirements. For example, when determining the 


'8 for example lunar or meteoritic material 


iia 


7 oh iv - 


a - 
a ee 
; oy : 
aley: ene nw iWesu ves ghey ea a ind iy he ed 
meet ~ : _—er 
eisviane @L by rb bt 240 Ag tay i. = 


Vi he 
noi taaty.« Peace ollie) ae Mbp buen: Ai aa 
‘gga sg? AAS) Oa) vain’ ee 


sed ar isem ¢en}as a ante add \a08 20 ile 


‘ yy ay ham ie 7¢ nt d2tep = 
Jf | Glenel? enw leakeggreant _™ 
je 3 Jivrna® Vitet aden AC a0 


. 7 
eeveaic 7 nie taey e437 16 


— 
nen = - ea> aif ct ,Wtaguate 307 


; <V9iw \amknaty BAS sane O0Ig 
‘ os Te3g : ‘ A beces in of 8 ao Lge 
nergzuwe Geyelsa be tamojye | 
_ neonecaitetni Ru woah setaris ees | jaeeee ass 26 a 
| i ] ad tn ,tevevrd , 3669 wal Bi 
ci 73 ; AS une ‘SR1gese a y SPAtaes en Ae sciamysitle 
5 


ree UCe)’ Pe be Lane vree gine 2t wt teres 
ine Jw Jaen ye ad (Gna elagtons SBR 
29 4 ee 

i (ew eh wol ‘pasyv ov a7 aes oy me 


19"S S65 «fen oor pid ileleta ve clhesats hy 
40%e  anaeFal Oem eo Slvone beara oases adkanlt 


pil * i os (2essResea @ '7H89 A pevnwed al 74 9 ae 
; a 7 


ad G2 ate 5 sao7 Oe ©. 619 La a ea 
va 7 


= ima ‘ame ute. eta 4 | ie 
- — a 


_ 


Le 


concentrations of the major heat-producing elements U, Th 
and K in rocks, K would have to analyzed by an alternative 
method. For these three elements, the increased speed with 
which they may now be determined using a sodium-iodide 
detector and the lower energy part of the spectrum, makes 
y-spectrometry an attractive alternative to the expensive 
(approximately $60 CDN/sample 1982 prices) XRF method of 
Galson et al. (1982). The former method suffers from the 
disadvantage that disequlibrium in the uranium chain may 
lead to spurious results, while XRF has the advantage that 
many others elements may be determined from the same pressed 
pellet with some extra analysis time. If uranium is 
determined by DNC or INAA in conjunction with natural 
y-spectrometry, this may give useful information regarding 
the state of disequilibrium in the ***U decay chain and 
hence the levels of the intermediate daughter products may 
be calculated. With ‘conventional’ XRF however the accuracy 
and precision of uranium and thorium analyses are poor (see 
Enoqurese te anGaai 6) - 

INAA, as has been shown, is a multielemental technique 
and extremely Sensitive for a large number of elements. Once 
again access to a reactor and gamma-ray spectrometry system 
usinecessary. For optimum Vh sensitivity, decay periods of 
two weeks or so are required making the method more 
expensive than say DNC. Because such small Samples are 
generally used for INAA care must be exercised in sample 


preparation to avoid sample inhomogeneity. 
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When comparing various analytical techniques for 
determining U, Th and K in granitic rocks Stuckless et al. 
(1977) concluded that y-spectrometry was adequate for the 
concentrations normally present in such rock types, but as 
disequlibrium in the ?**U decay chain was occasionally a 
problem, they suggest that the method be combined with DNC 
for uranium analysis. Hart et al. (1980) found when 
comparing INAA, XRF and mass Spectrometry isotope dilution 
(MSID) for U and Th analysis that at concentrations above 10 
ppm U and Th, INAA and XRF are alternatives to the expensive 
and time comsumming MSID technique. As a wide variety of 
geological samples will have lower concentrations than 10 
ppm U they state that INAA is to be preferred over the less 
precise XRF. 

The uranium results of the various techniques employed 
in this#work are listed in table 7.1 and ane shown 
GEapnically sin £19Gures Viet ow.5. OF tEhevmelnods suri tized 
XRF and FTA are the least precise'’. The agreement between 
the INAA short- and long-irradiations, and the 
gamma-spectrometry and DNC results is most satisfactory. 
Figures 7.6 and 7.7 show the comparison of thorium analysis 
by low energy gamma-spectrometry and XRF with INAA 
determination of thorium. 

The choice of which method to use for uranium analysis 
therefore depends ultimately upon the user requirements and 
it is hoped therefore that this part of the thesis on 


‘2 in the case of FTA the imprecision is largely due to 
hetrogeneities in the comparator standards. 
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Table 7.1 Results of various methods of uranium analysis 
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applied to Exshaw shale samples 


(2)* gamma-spectrometry with no potassium correction 
(3)* gamma-spectrometry with potassium measurement and correction 
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instrumental analytical methods of uranium, thorium and 


potassium analysis will be of use to the prospective 


analyst. 
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Figure 7.1 Linear plot of uranium content obtained by DNC 


versus XRF 
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Figure 7.2 Linear plot of uranium content obtained by DNC 


versus low energy gamma-spectrometry NaI analysis 
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Figure 7.3 Linear plot of uranium content obtained by DNC 
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Figure 7.4 Linear plot of uranium content obtained by DNC 


versus INAA 4 minute irradiation 
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Figure 7.5 Linear plot of uranium content obtained by low 
energy gamma-Spectrometry with and without correction for 


potassium 
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Figure 776 Linear plot of thorium content «obtained bya2 hour 


INAA versus XRF 
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Figure 7.7 Linear plot of thorium content obtained by low 
energy gamma-spectrometry Nal, corrected for potassium 


versus XRF 
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Geochemistry of the Exshaw Shale. 
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8. The Exshaw Shale 
The Exshaw Formation, although originally thought to have 
been Devonian in age (Warren, 1937), was dated by Folinsbee 
and Baadsgaarde 1950) ate267e2e139My.. making it 
MissisSippian (Lower Carboniferous) in age. 

The Formation is marine in nature and consists 
essentially sof twomunits in Alberta =van upper si bry 
limestone and a lower black shale (Havard, ibid) overlying 
the Wabamum Formation on the plains and the Palliser 
Formation in the Rocky Mountains and Foothills. On a 
regional scale, the Exshaw Formation is an extension of a 
continent wide Devonian - Mississippian black shale sequence 
and as early as 1948 was correlated with the uraniferous 
Chattanooga Shale of the Mississippi Valley and 
Mid-Continent regions of the United States (Fox, 1948). 

Ettensohn and Barron (1981) proposed a depositional 
model to explain the extent and conditions responsible for 
this continent wide black shale sequence. They suggest that 
the Acadian Mountains, produced by continent-continent 
Collision, formed a north-south trending barrier vacross the 
equator on the east side of the Old Red Sandstone Continent. 
This mountain range acted aS a barrier to moisture-laden 
trade winds, producing rain shadow conditions west of the 
mountains were two epicontinental seas, Separated 
intermittently by the Trans-Continental Arch, were located: 
sea result of being Imvthis rainshadow area, “the reduced 


clastic influx and restricted circulation within the basin, 
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led to the deposition of black, euxinic, organic-rich muds. 


8.1 Mineralogy 

In hand specimen the shales are fine-grained, fissile 
and vary in colour from black to dark grey, depending upon 
their organic and carbonate content. Examination of a number 
of back-packed, randomly orientated bulk samples was carried 
out by X-ray diffraction (XRD) analysis to determine sample 
mineralogy and to see if there was any obvious mineralogical 
control on the U, Th and K contents of the samples. To this 
end a Philips wide-angle goniometer-diffractometer and strip 
chart recorder were used. Samples were selected from their 
USSTR and kKycontents anderunsgover van26 trangerof 2-604 at 
constant instrument settings. This allows the direct 
comparison of peak intensities between samples so that 
relative concentrations of the different mineral phases may 
be obtained. Cu Ka radiation was used for all the 2-60° 
runs. However, due to a variable 'spike' at the lower 26 
angles in these runs, the majority of the samples were 
re-run over 2-20° 28 on a second Philips system using Co Ka 
radiation, 

The 'd'-spacing for each major reflection was 
calculated and from them the minerals identified. 
Reflections due to quartz, calcite and dolomite dominated 
the traces, although those of orthoclase and muscovite were 
also present in some samples. Although not eaSily identified 


(as no clay seperate was taken for analysis) illite was 
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found to be the dominant clay mineral group present in the 
samples studied. This is in agreement with Havard ibid who 
carried out more extensive XRD analyses on Exshaw Shale 
Samples from southern Alberta, and various other workers who 
found illite to be the dominant clay group present in 
Palaeozoic black shales (e.g. Bloxam, 1964; Leventhal and 
Hosterman 1982). 

No obvious mineralogical @eontroel waseseensforsU and Th. 
However, aS might be expected, K enrichment seems to be 
related to the orthoclase and illite or muscovite content of 
the shales. 

The relative amounts of dolomite/calcite versus illite 
may possibly be used as a broad guide for differentiating 
the Exshaw Formation into its two divisions: the upper silty 


limestone or dolomite and the lower black shale unit. 


8.2 Geochemistry 

Ahrens (1954) noted that in most instances the major, 
minor and trace elements of samples taken from the same 
population displayed a log-normal distribution. Histograms 
of the results of this work were found to either follow a 
lognormal distribution or to closely approach @t within the 
limits of the number of samples analyzed. This is important 
aS Many parametric statistical tests commonly employed (e.g. 
correlation coefficients, cluster analysis etc.) assume a 
normal distribution. Consequently, before any statistical 


analysis of the data could be attempted using these tests, 
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log transformation of the data had to be completed. This, 
and subsquent statistical analysis of the data on the Exshaw 
shale, was carried out using the Michigan Interactive Data 
Analysis System - MIDAS (Fox and Guire, 1976) run on the 
University of Alberta Amdahl 470/6 computer. 

In an attempt to determine and aid in explaining the 
relationships between the various elements analyzed Pearson 
product-moment correlation coefficients were determined for 
all pairs of variables. This coefficient is a measure of the 
strength of the linear relationship between any two 
variables in queston, that is, how well the variation in one 
element can be explained by that of another. 

The samples were first arbitarily divided into 'black' 
and 'calcareous' shales based on their CaO content. Those 
samples with more than 5.0% CaO being classified as 
‘calcareous'. This wasS carried out so that any subtle 
correlations would not 'cancel' each other. 

The median values for each element determined in the 
course of this study are listed in table 8.1 together with 
"'average' shale values from Turekian and Wedepohl (1960) and 
metal-rich black shales taken from Vine and Tourlelot 


(1970). 


8.2.1 Major elements 
Suen. 
The median silica content for all the Exshaw samples is 


48.5% SiO, while the calcareous and black shales give 36.9% 


vai 7 


(oie? (hespipentt ae ‘cache 


’ ' 


vadesd att im anae ait Ye pivytion sebtowtiade: Saat = 


nr neh rem 


gall ev east) re@Zns if ots enity wor steel Ladle 


hs er 


atid owhloe Bim wot) G20 © oocgeyhs 


[- 
7 SATs avers sitetseA Ai sdia ja efi ewialy : 


5 ‘ (Av — » wy Lease | ae) 3a AG at a 7 
>. Le 
newcod egidencizelai 7 
_ 7 
i 4 > Apse T79CR MG 4.416% 08 ma dwar : 


ie y iy azheq bie : 


. : _ 
s Gas ® :2 2 Aspnerge _ 
f _— 
4 2; i - Vol a 48 ea Loe | aeldsiseyv. a 
__ 
- 
aigas ed nam Jnemel® . 
- _ 
sn | eh lw etbunsiitu. 399i. een ae am. oT 
’ ' : ; ¢ f vos eotan’ a 
’ (can etew Avie aniquaé 
- 
: = 
a » | A t » uns g ‘sg he e soles’ 
7 
¢ Slute enotseie7309 
\ : f ~w.8 : , seutav @albea sft a 
seal eodhanos 7 eae se ylege eine Yo eahude 
hai ' jehew bin Aaide aqisev ataca pau 


37¢@{,@A6 isa ié steel im, 


el be (yao T the a als nid 
oe ants wl pe 


128 


table O21) Comparison of average shale values with those of 


the Exshaw Shale determined in this study 


Element Average Black Calcareous Average* Average ft 
Exshaw Shale Shale Shale Shale black shale 
Sid, 48.5 58.6 39 .6 5Sielt = 
Al,0, 7.20 9.65 5.45 NGA No 
Fe,0, 3.88 4.46 3.43 675 2.86 
MgO Io eAl 0.84 He git 2.49 1.16 
Cao 4.92 1.69 13.0 3.09 2.10 
NaO 0.30 0.34 O25 es) 0.94 
K,0 2.67 SoS 2.16 £570 2.41 
TiO, 0. 36 0.42 OF 0.77 0.33 
As 21.0 27.8 1633 ig} == 
Ba 325 355 299 580 300 
Gil 470 565 394 180 SSS 
Co 20.6 24.1 1759 19 10 
(ere B58) 66.0 47.0 90 100 
Cs 4.08 oes} SAS 5 === 
Cu 72 82 62 45 70 
Ga 1253 WHS eal 19 20 
Hf DB 28) 19 2.8 ——— 
Mo 27.6 21.6 Ban0 2.6 10 
Ni 107 123 92 68 50 
Pb 26 30 22 20 20 
Rb 78 101 61 140 a 
Sb 3.56 4.36 2eo8 ss = 
Se Jealal IYO) 38) 7.74 13 10 
Sr 222 157 303 300 200 
Ta 0.65 0.76 OS 0.8 Sss5 
Th oZ 10.7 8.2 WZ a 
U 15.6 1529 N52 Biv) — 
Vv 227 267 TOM 130 150 
an 102 108 M5) 95 30 


* Average shale values from Masan (1966) after Turekian and Wedepohl (1961) . 
t Average black shale values form Vine and Tourlelot (1970). 


(Major elements in weight percent, trace elements in parts per million, ppm.) 
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and 39867) "respectively. The black shales ishow "a tstrong 
negative correlation with Al(r=-0.86) and Ti (r=-0.86) which 
1S contrary to Havard's findings. He notes however, that 
much of the silica present was biogenetically precipitated 
as sponge spicules. Therefore, it would be surprising to 
find that Si would covary with Al and Ti except by chance, 
when much of the silica present has a completely different 
Orlgium to ‘that “of “the clay. The “silica ‘content ef <the 
calcareous shales shows a strong negative correlation with 
Ca (r=-0.82) which as expected shows that conditions 
conducive to limestone formation are very different to those 
Suitable for black shale formation. 

Pron. 
All the samples gave a median Fe203 content of 3.88%, while 
the black shales give 4.46%. The Exshaw black shale value is 
greater than the average black shale concentration listed by 
Vine and Tourlelot (1970). Except for a moderate correlation 
with AS (r=0775) forsall’ the samples (0076"tor the black 
shales) iron was not found “to “covary “Significantly ewith any 
other element. Iron pyrites, macroscopically present in many 
of the samples (particularly the black shales), obviously 
accounts for much of the iron present. Arsenic, which is 
concentrated at more than twice the average shale value in 
the black shales, is presumably incorporated largely within 
the pyrite structure. The median Fe values obtained are 
greater than those reported by Havard. This may in part be 


explained by the fact that a number of his samples were from 
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Surface exposures and consequently it may have been 
difficult to obtain fresh, unweathered samples. 

Aluminium. 
Good agreement is found between the median Al,0; values of 
this work and that determined by Havard. For all the samples 
Al shows strong positive correlations with Rb (r=0.90, 0.98 
EOmache black shales alone), Ga (r=0066)', Ti (r=0,62), Sc 
VEe=C eee Ka CrH0079) 7 Cr (CeO. 76)" Th GeO oaimand Ts 
(r=0.76) (see figures 8.1 and 8.2) and to a lesser extent 
with Hf, As, Ce, Cs and V. The majority of these are 
hydrolysate elements (see figure 8.3) and are either 
structurally combined within the clay minerals and/or are 
adsorbed by them. The correlations with the soluble cations 
Rb, K and Cs, suggests that much of the Al also occurs 
within potassium feldspars, both detrital and authigenic in 
origin. The poor correlation between Al and the REE 
suggests, together with other information given later, that 
the REE's are associated with the clay minerals and that the 
low REE content of feldspars serves to 'dilute' the 
correlation with the clays. The negative correlation between 
Al and Ca (r=-0.63) in the calcareous shales supports the 
Silica-calcium carbonate conclusion that limestone and black 
Shale formation are incompatanle. 

Calcium and Magnesium. 
The median values of CaO and MgO differ markedly from the 
black to calcareous shales, namely 1.69% to 13.04% CaO, and 


0.84% to 1.71% MgO, respectively. As shown from some of the 
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Figure 8.3 Diagramatic representation of ionic charge versus 


ionic radius of elements in the hydro-geological sedimentary 


environment 
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XRD analyses a number of samples contain calcite or dolomite 
and clearly these are the major sites for the two elements. 
There 1S practically a random relationship between Ca and Mg 
in the calcareous shales in agreement with the general 
observation that calcite and dolomite do not often occur 
together. 

Potassium. 
The median K;0O content varies from 3.32% for the black 
shales, which is close to the 3.24% of the 'average' shale 
and in good agreement with Havard's 3.29% for all his 
samples, to 2.16% K20 for the calcareous shales. The 
strongest correlation observed is with Rb (r=0.79 for all 
Samples, see figure 8.4) which is not unexpected on account 
OfVtheiresimilar™tonie radii and identical charge. Brom the 
correlations with many of the other hydrolysates and soluble 
cations potassium, like aluminium, is thought to be located 
in both the clay minerals (dominantly illite) and in 
detrital and authigenic potassium feldspars. 

Sodium. 
Of particular note is the very low concentration of Na2zO 
(0.32% for all samples) compared to the average shale value 
Of sie30%, Silrahtlysmore thane four times lower echnan ene 
average shale this emphasises the dominance of illite as the 
major clay group present (as opposed to the more sodic 
smectites), and potassic detrital and authigenic feldspars. 


The little sodium that is present probably does occur with 


thesclay mineralsi. 
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Figure 8.4 Log-log plot of K2O versus Rb for Exshaw Shale 
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Titanium. 
The median TiO, value for all the samples is 0.36%. Although 
in excellent agreement with Havard's average value of 0.35%, 
it 1s much lower than the average shale value of 0.76%. The 
Surong scorrelation with Al -andiK,rtworot the major clay 
forming elements, together with the hydrolysate elements Sc, 
Cr, Ga, Th, Ta and Hf show titanium to be largely associated 
with the clay fraction and not, as is often the case, only 
with detrital resistates such as ilmenite and sphene. Hirst 
(1962) noted from studying the geochemistry of modern 
sediments in the Gulf of Paria that the Al,0;/TiO, ratio of 
the sediments was fairly constant and that it was unlikely 
that much Ti waS precipitated out of solution. He suggested 
that it was more likely that the Ti was associated with the 
clay fraction during weathering and transportation. As the 
Exshaw shale sedimentation rate is thought to be have been 
very low, and there is a lack of coarse clastic material 
present, the low Ti content may therefore be a consequence 


of the paucity of detrital material. 


8.2.2 Rare Earth Element Geochemistry 

As noted earlier, in the process of analyzing samples 
for uranium, thorium and potassium by INAA a number of other 
elements may be determined simultaneously. Depending upon 
the rock type and matrix, seven to nine of the fourteen rare 
earth elements (REE) may routinely be quantified by the 


procedure described. Due to their similar ionic radii and 
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preferred trivalent oxidation state the REE generally behave 
coherently*°®, and consequently they have been the basis of 
many geochemical studies. Traditionally, they have been 
utilized in hardrock research, though over the past two 
decades their application to sedimentary problems has been 
on the increase. 

Overview of REE in sedimentary rocks 
From the early work of Haskin and Gehl(1962), Haskin and 
Frey ( 1966))"Haskin®et al. “%1966)) and °Ronev’ etal: 
(1967,1974) a number of features became apparent regarding 
the REE contents and patterns of sedimentary rocks. Of 
Particular note was the similarity of REE patterns (when 
normalized to chondrite REE values) between shales of 
various ages and locations (albeit that the relative 
abundances differed). This helped compound the idea that 
during shale formation exposed upper crustal material was 
both representively sampled and mixed in an efficient 
manner. 

In addition it was found that the total REE contents in 
elays®> Sands > *carbonates, ana™ that the REE ares generally 
adsorbed by clay minerals and transferred almost 
GUantitatively intorclastic Sediments. Relative to 
chondrites, sediments are enriched in all the REE, with an 
increasing enrichment of the light REE (La~-Sm). 

Due to the similar REE patterns of sedimentary rocks 
Haskin and Frey (1966) published the results of the analysis 


2°Ry and Ce are exceptions, Eu** may be reduced to the 
State ana Cem oxidized to the = state; 
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of a composite of forty shale samples, called the North 
American Shale Composite (NASC) with which most sediments 
and sedimentary rocks have been compared (see for example 
Haskin eb ealva(966), and Dypvik and Bruntelt (1976, 1979) ). 

As REE contents and patterns of clastic sediments were 
thought to reflect the parent material of the sediment, 
Ronov et al. (1967) suggested that there may be a major 
change in the patterns and absolute abundance of the REE at 
the Archaean-Protoerozic boundary, when major additions of 
felsic material to a predominantly mafic continental crust 
were thought to have occurred. Their few analyses seemed to 
Support this idea. Work by Nance and Taylor (1976), and 
McLennan et al. (1980) showed conclusively that such a 
pattern change does infact occur. Analysis of essentially 
unaltered Archaean (world-wide distribution) and 
post-Archaean sedimentary rocks (from Australia) showed the 
former to have no negative Eu anomaly when chondrite 
normalized compared to the post-Archaean Australian 
sedimentary rock (PAAS) and the NASC. In addition analysis 
of the various shales of different ages that made up the 


PAAS showed there to be an increase in the absolute 


abundance of the REE with time (figure 8.5). 
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Figure 8.5 REE normalized patterns for Australian shales of 


various ages 
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Figure 8.6 REE normalized patterns for various sedimentary 
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Results of this work 

The results for each Exshaw Shale sample analyzed are 
given in Appendix 1. The median and extreme values, 
normalized to chondritic values (Nakamura,1974) are plotted 
in figure 8.7, together with the North American Shale 
Composite, Post Archaean Australian Shale and the Alun Shale 
of Norway (Dypvik and Brunfelt, 1976). The Alun Shale is 
lower Palaeozoic in age, slightly uraniferous 
(Bjgrlykke,1974), is thought to have been deposited ina 
shallow epicontinental sea, and the samples analyzed: 

"represent typical black shales high in 

organic carbon and sulphides, deposited under 

extremely Slow sedimentation rates, about 

imm/1000 years." Dypvik and Brunfelt (1976) 

Be 2a 
Clearly then the conditions under which both the Alun and 
Exshaw Shales formed were very similar, and this may allow 
analogies to be drawn between them. In addition to the 
conditions under which they formed, their REE patterns are 
quite Similar even to the extent that they both have 
negative Ce anomalies (see figure 8.7). 

To compare the mean Exshaw results with other shales 
the median values have been normalized to the PAAS and the 
results plotted in figure ¢ 63 Erom this) dtaqnuamn, 1t can be 
seen that the total REE of the Exshaw Shales are less than 
those of both the PAAS and NASC. This is in part due to the 
older age of the Exshaw Frm. compared to many of the samples 
making up the PAAS (see figure 8.5) but more importantly it 


Ps) thought to be due to dilution of the shales by carbonates 
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Figure 8.7 REE normalized patterns of the Exshaw Shale and 


comparatory shales 
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Figure 8.8 REE pattern of the mean Exshaw Shale normalized 


to the post-Archaean Australian shale 
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(see figure 8.6) and possibly by biogenetic silica (see 
later) both of which are known to be depleted in the REE 
relative to clays (Nance and Taylor (1976) and Piper 
(1974)). A similar effect was observed in some of the 
samples from the Oslo epicontinental sequence, where it was 
Suggested that depletion in the total REE contents of two 
samples, 56 and 57 was: 

"probably the result of dilution by relatively 

Marge ameountsSeot calcite in aso and <quastz in 

5) Deo ho “Dy pV MepandeBuuntele 197-6 9:. 
The smaller Exshaw Eu anomaly suggests that there is less 
feldspar in the Exshaw than in the PAAS, as feldspar 
concentrates Eu during magmatic processes. From the PAAS 
normalized plot Sm and the heavy REE (HREE) Dy, Tb, Yb and 
Lu are quite smooth (although depleted relative to the PAAS 
by 10 - 20%) while the light REE (LREE) La and Ce are 
Significantly depleted relative to the PAAS. Of particular 
note is the negative Ce anomaly. 

What is the cause of this negative anomaly in the 
Exshaw Shale? Ronov et al.(1967) noted that Ce** may be 
oxidized to the Ce** in zones of weathering and 
sedimentation, and that in this oxidized state its behaviour 
is more similar to the HREE than the LREE. In addition, of 
all the REE only Ce may be oxidized in seawater to the 


"highly insoluble quadrivalent state” (Piper 
1974, p.286) 


Piper (1974) found that Ce was significantly depleted an 
seawater as well as in foraminifera (after Spirn, 1965) and 


diatoms, these organisms inheriting their REE patterns from 
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seawater. Nance and Taylor (1976) analyzed and plotted the 
REE pattern of a Carboniferous limestone which displayed a 
negative Ce anomaly and Ce/La ratio of 1.68 (cf. NASC and 
PAAS ratios of 2.28 and 2h 10 sespectively)¢ sfheyedia not 
however find a Ce anomaly in the PAAS and explained it as 
follows: 

"The absence of cerium anomalies observed in 

this study indicates that equilibrium during 

deposition with seawater, which has a 

pronounced negative cerium anomally, does not 

takesplace” (Nance and! Taylor; 976)- p. 1043) = 
A number of factors appear therefore to have contributed to 
the negative Ce anomaly of the Exshaw Shale. Dilution of the 
black shale component by the calcareous shales appears in 
part to have been responsible, as is shown by the chondrite 
normalized plots of the two shale components in figure 8.7. 
The extremely slow sedimentation rate of the Exshaw Shale 
would also have offered greater opportunity for equilibrium 
to have been approached, if not gained between sediment and 
Seawater, than would more rapid deposition. This idea is 
supported by the extreme similarity of the REE contents of 
the Exshaw black shales and the sedimentary rocks from the 
Qslo region of Norway (Dypvik and Brunfelt, 1976). Finally, 
as mentioned earlier, much of the silica present in the 
Exshaw Shale is of a biogenetic origin (Havard, 1967). 
Presumably the silica of sponge spicules precipated from 
seawater would also have a negative Ce anomaly similar to 
that found for diatoms (also siliceous) and foraminifera.?? 


21 'anomally' spelt as reported in publication. 
22 The author knows of no REE analyses of sponge spicules to 
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Towards che end’ Of this study the publication by Tlig and 
Steinberg (1982) regarding recent sediments of the Indian 
Ocean added weight to this idea. One of their findings was: 
"The negative Ce anomaly of the coarsest 
fraction is associated with the biogenic 
components (foraminifera, diatoms) which 
inheret their REE distribution from seawater." 
(Tlig and Steinberg (1982) p.331) 
However, it 1S not possible to explain the observed anomaly 
soley by dilution of an 'average' shale, say the NASC or 
PAAS, by sponge spicules. It is therefore concluded that the 


anomaly 1s the combined result of dilution by carbonates and 


Lie conmteribution Of biogeneeic: fractions. 


REE and other element correlations 

As one would expect, all the REE are strongly 
correlated with each other. For example the correlation 
between La and Ce is 0.96, between Dy and Sm 0.98 and 
between Tb and Sm 0.93 (see figures 8.9, 8510, 8.11). Most 
of the remaining strong REE correlations are with 
hydrosylate elements, in particular Hf, Ta, Ga and Th which 
agrees with the hypothesis that the REE are largely 
associated with the clay minerals. The correlation between 
Th and the REE (in particular La) is of particular interest 
and is the subject of a publication by McLennan and Taylor 
(1980). They note that the La/Th ratio of mafic igneous 
rocks WG greater than thateot felsic ei gqneoussrocks: 


Statistically they were able to show that the La/Th ratio 


22(cont'’d)support this idea. 
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Figure 8.9 Log-log plot of La versus by for Exshaw Shale 
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Figure 8.10 Log-log plot of Dy versus Sm for Exshaw Shale 
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Figure 8.11 Log-log plot of Dy versus Eu for Exshaw Shale 
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for the Australian Archaean Shale (AAS) is greater than that 
of» the PAAS and suggest that this ts additional evidence for 
a more mafic upper crust in Archaean times than at present. 
some 75% of the PAAS analyses lie between the La/Th ratio of 
2 and 4. Results of this work show greater variability, 
which possibly suggests that during Palaeozoic times the 
exposed Canadian Shield was quite variable in composition. 
Of interest are the average La, Th, U, Rb and K20 
values for the Archaean and post-Archaean crust calculated 
by McLennan and Taylor (1980). These values, together with 
the median values for all the Exshaw Shale samples and the 
black shale samples are listed in table 8.2. Except for the 
anomalous uranium values, the agreement between the black 
Exshaw Shale samples and both the PAAS and the average 
values for the Canadian Shield (Shaw et al. 1967, 1976) are 
excellent. This agreement suggests that black shales may be 
used to estimate the average composition of the Source area. 
On a regional scale they may possibly therefore, be used to 
evaluate the average composition of the continental crust. 
As suggested by McLennan and Taylor (1980, 1980b) it would 
then be possible to demonstrate whether the composition of 
the continental crust has been changing with time. Once 
again dilution of the shales by the carbonate component may 


account for the lower median Exshaw Shale values. 
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Table 8.2 Comparison of some elemental concentrations in the 


post-Archaean upper continental crust with the Exshaw Shale 
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Table 8.2 Comparison of some elemental concentrations in the 


post-Archaean upper continental crust with the Exshaw Shale 


K50 (3) Sie 35350 Sigil Ze O" S232 


(after McLennan and Taylor,1980) 


a. The post-Archean upper continental crust (McLennan et al.,1980) 

b. ‘The post-Archean upper continental crust (Taylor and Mclennan, 1980) 
c. ‘The average exposed Canadian Shield (Shaw et al., 1967, 1976) 

d. The median of all Exshaw Shale samples analyzed (this study) 

e. The median of all black Exshaw Shale samples analyzed (this study) 
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8.2.3 Uranium and Thorium geochemisty 

Uranium 
Maximum and minimum U values for the Exshaw Shale samples 
rangenfiromi92.2eppmrTeo 3.4 ppm? swith atmedtan valuerotel5.6 
#o2,2 ppm. (This*value is in agreement withethe 136+ 8 ppm 
mean and standard deviation reported by Campbell (1980) for 
the Exshaw Shale from a frequency vs. log uranium plot. The 
15.6 ppm is a factor of four greater than the ‘average' 
shale value of 3.7 + 0.5 ppm (Adams and Weaver, 1958) anda 
factor of five greater than the average post-Archaean shale 
at 3.1 + 0.2 ppm (McLennan and Taylor, 1980b). The level is 
however less than the Chattanooga Shale median uranium value 
of 20 ppm (Landis, 1962). In a given core large variations 
in the uranium concentrations were found, for example, one 
core showed a minimum and maximum concentration of 12.9 and 
91.4 ppm U. Variations, similar in magnitude to this have 
been observed in Devonian black shale cores by Leventhal 
(1928 ).. 

The median values for the black and calcareous shales 
at 15.9 and 15.2 ppm U respectively show that there appears 
to have been no obvious mineralogical control in either 
Group whieh Controlled sthe enrichment. 

A number of workers (Bloxam, 1964, Mo et al. 1973, 
Kelepertsis, 1981, “heventhaly 91978 sand Leventhal and 
Hosterman, 1982) found there to be a strong correletion 
between uranium and the organic carbon content of black and 


calcareous shales. This suggests that uranium was 
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precipitated out Of solution oy the reduction of the Soluble 
Us lon to’ the insoluble U** ton which was then e:ther 
adsorbed by, or complexed with, the organic material 
Present. The tabulated U, Th, K2O and organic @ correlation 
coefficients for the Exshaw Shale samples are given in table 
8.3. Uranium is most strongly correlated with organic 
carbon, particularly in the black shales. This suggests that 
uranium waS supplied in solution (as it has only wea 
correlations with most of the remaining thirty or so 
elements determined) from a low lying source area and 
precipitated out of solution in a reducing environment under 
conditions suitable for the preservation of organic matter, 
and was able to attain the concentrations found due to the 
very low sedimentation rate. 

The only other element which uranium is found to 
Correlate Significantly with is mickel (0.70). This finding 
is in partial agreement with the sympathetic variation of U 
with Ni and Zn observed in the Exshaw Shale by Campbell 
(1980). | 

Thorium 
In the Exshaw samples analyzed thorium concentrations ranged 
from i123 toe43.8 ppm, with a median value ote 32274 sic 9 pom 
(N=42) 1 sigma). This is Slightly lower than the average 
shale value of 12 + 1 ppm Th reported by Adams and Weaver 
(7959), or the 12.3 ppm Th average concentration reported by 
Bloxam (1964) for Carboniferous black shales from South 


Wales and is less than the post-Archaean sedimentary rock 
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Table 8.3 Correlation coefficients for U, Th, 420 and 


organic C for the Exshaw Shale 


All Samples. 


Black Shales. 


U 
Cal 
a careous 
Shales. 
K,0 


(Sishe(sig 


average value of 14.6 + 1.2 ppm (2 sigma) (McLennan and 
Taylor,-1980b). The median black Bxshaw* Shale concentration 
of 10.7 + 2.1 (N=19) is however in better agreement with the 
average shale values. The respective calcareous and black 
Shale wmedian thoriumavalues of 8542 2. lc 3eppme(N=23) sand 10.7 
+ 2.1 ppm (N=19) show there to be little or no significant 
difference between the Exshaw calcareous and black shale 
thorium concentrations. Kelepertsis (1981) found the average 
Th content of intercalated Carboniferous calcareous and 
black shales in Anglesey, North Wales to be 16 + 2 ppm (N=5) 
and 11.5 + 1.3 ppm (N=4)respectively.?°* However, due to the 
few samples analyzed little significance may be attached to 
his findings. 

From the earlier disscussions regarding LREE-Th 
relationships and the correlation coefficients between Th 
ana sk.O-eShowm.im table G23) Sht.iS apparent sihat emuch:= 15 not 
all, of the Th in the Exshaw Shale samples iS associated 
(adsorbed) with the clay minerals. The strong correlations 
with Al (0.79), Hf (0.90) and Ga (0.82) supports this 
hypothesis. Bloxam (1964) found that Th and K exhibited some 
degree of covariance and suggested, as did Pliler and Adams 
(1962) and Kelepertsis (1981) that Th was most likely 
associated with the clay fraction (usually illite) of the 


shales they studied. 


23Kelepertsis (1981) actually reports 19 and 12 ppm Th © 
respectively for the calcareous and black Shales but failed 
to log-transform his data before analysis. 
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AS This extremely insoluble im near neutral ‘waters 
(Kelepertsis, 1981) the slightly lower than average Th 
content of the Exshaw samples agrees with the hypothesis 
that during deposition of the sediments, detrital influx was 
low as “the Th was most Zikelly incorporated in detrital 


clays. 


8.2.4 Geochemistry of remaining elements 
Alkali metals and Alkaline Earths. 

Particularly in the black shales, cesium is found to 
correlate with Rb and K while Rb is strongly correlated with 
Ae?) Ga (0590) Ba heeSewand Cn. sardine eatkeomroune 
to correlate strongly with these elements. This suggests 
that Cs, Rb and Ba are largely associated with the clay 
minerals. Strontium however, is clearly associated with the 
calcareous shales. This association was also noted by Havard 
ibidvand more generaily in Shales by Vinewands Tourlelor 
(1970). 

The Base Metals - Cu, Ni, Pb, Zn and Co 

The transitional metals and Pb are all enriched in the 
Exshaw black shales compared to the average shale values of 
Turekian and Wedepohl(1961). With the exception of Pb, no 
Significant Correlations are observed between these metals 
and any other element. When enriched these metals are often 
found to be strongly correlated with organic carbon, with 
which they often form organo-metallic complexes (Leventhal 


and Hosterman 1982). Lead is correlated with the LREE, Th 
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and many of the other hydrosylate elements suggesting that 
it is associated with the clays. The correlation with Th 
(0.64) agrees with the findings of Weiss and Amstutz (1966), 
who noted’ that. Pb?- “could interchange with Th** in clays. 
Their experimental studies showed that between 2 - 6 weight 


percent Pb may be contained within the illite structure. 
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9. Conclusions 


9.1 Methodology 

Of the various methods of uranium analysis used the 
delayed neutron counting method is to be favoured if only 
uranium is to be quantified and sufficient sample is 
available. Where only small quantities of sample are 
available, fission track analysis or INAA are sensitive 
alternatives to DNC. 

When DNC is used in conjunction with the gamma-ray 
spectrometry of naturally occurring radioisotopes cases of 
5% disequlibrium or greater can be identified. Where, in 
addition to uranium, multielement information is required, 
INAA 1S to be prefered. Conventional XRF analysis for the 
quantification of U and Th although reasonably accurate 
lacks precision at typical rock concentrations. 

From the analyses of the great variety of rock types 
analyzed (standards and samples) INAA has shown to be an 
accurate, and in most inStances precise, method for use in 


the earth sciences. 


9.2 Geochemistry 

The anomalous radioactivity of the Exshaw Shale is 
largely due to uranium enrichment (up to 92 ppm), though 
potassium levels are in certain instances high (maximum of 
almost 8% K20 by weight). The median thorium value for all 


the Exshaw Shale samples is less than the 'average' shale 
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value. 

Uranium appears to have been transported in solution, 
and upon reduction to the insoluble tetravalent state, was 
either adsorbed or complexed with the organic matter 
present. In most instances the concentrations of elements 
from the Exshaw Shale agree well with the values from 
Turekian and Wedepohl (1961). 

Exshaw Shale rare earth element patterns show a 
negative Ce anomaly in many of the samples analyzed; the 
anomaly being the result of dilution of the black shale 
component by calcareous matter and by biogenetic silica. 

Agreement between the median elemental concentrations 
of the black Exshaw Shales and the average values of the 
Canadian shield (Shaw et al., 1967, 1976) suggests that 
during black shale formation variations in the source area 
geochemistry are throughly mixed. Consequently, median 
values, determined from many Shale samples, or froma 
composite shale sample, may be used to estimate the average 
composition of the exposed source area, and if over a large 


enough area, the average composition of the exposed 


continental crust. 
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Appendix 1 : Chemical Analyses. 


the winstrumental neutron activation analysis and X-ray 
fluorsence results for the forty-two Exshaw Shale samples 
ene sreported in tnismappendix, Sieicarseu, Ni, sPbrandecn 
webe analyzed by ARE, organic carbon with a BECO high 
frequency induction furnace in a manner Similar to that 
outlined by Boyce and Bode (1972), and the remaining 
elements via INAA. 

The INAA resultS are reported with + one standard 
deviation uncertainty and total Fe is reported as Fe203. 
NA - element not analyzed for 


ND - element not detected during analysis. 
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Lu 1.17 + 0.03 0.37 + 0.04 0.22 + 0.02 
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Dy 5.21 + 0.12 3.04 + 0.08 4,38 + 0.10 
Yb Sa OH 2.34 + 0.29 2.45 + 0.21 
Lu 0.37 + 0.03 0.30 + 0.02 0.28 + 0.03 
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Feo0, 2.30 + 0.15 8.76 + 0.25 4.51 + 0.19 
MgO 0.14 + 0.08 0.72 + 0.19 2.42 + 0.15 
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Eu 0.91 + 0.06 5.93 + 0.15 1.58 + 0.10 
Tb 0.19 + 0.06 3.91 + 0.10 1.10 + 0.06 
Dy 2.43 + 0.06 21.6 + 0.3 5.24 + 0.11 
Yb 1.48 + 0.17 4.69 + 0.32 2.24 + 0,22 
Lu 0.20 + 0.02 0.54 + 0.17 0.38 + 0.15 
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$10, NA 54.72 67.92 
A1,05 10.8 + 0.2 15.1 + 0.2 4.49 + 0.07 
Feo03 4.64 + 0.64 4.34 + 0.22 3.23 + 0.16 
MgO 1.24 + 0.18 1.33 + 0.23 0.74 + 0.10 
CaO 9.33 + 0.10 1.26 + 0.12 2.06 + 0.11 
Na 0 0.209 + 0.004 0.754 + 0.013 0.375 + 0.007 
K0 4.48 + 0.30 3.79 + 0.18 1.54 + 0.18 
Tid, 0.44 + 0.03 0.63 + 0.06 0.19 + 0.02 
C org. Se 2as 7.4 

As 17.6 + 0.4 29.1 + 0.5 19.8 + 0.6 
Ba 358 + 120 464 + 44 108 + 12 
C1 325 + 20 2040 + 70 955 + 37 
Co 2752 49025 21s te0ne 2204 43053 
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Cs 7.97 + 0.59 8.67 + 0.19 2.53 + 0.12 
Cu NA 47 131 
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Ni NA 106 82 

Pb NA 25 21 
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Sc 11:3°+ 0:1 16.8 + 0.2 4.92 + 0.16 
Sr 535 + 57 85 + 29 56 + 20 
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Th 10.2 + 0.9 12.3 + 0.6 5.05 + 0.93 
U 14.8 + 0.2 5.72 + 0.15 24.5 + 0.3 
V 185 + 5 499 + 9 138 + 4 
Zn NA 143 15 

La 23.8 + 1.4 PEP ails 12.9 + 0.5 
Ce ee) 55.5 + 0.5 23.4 + 0.4 
Sm 2.06 + 0.14 5.58 + 0.35 2.33 + 0.25 
Eu 0.48 + 0.04 1.29 + 0.08 0.42 + 0.05 
Tb ND 0.20 + 0.05 ND 

Dy 1.53 + 0.04 4.05 + 0.10 1.66 + 0.06 
Yb sees Wy) 2.37 + 0.22 ND 

Lu 0.65 + 0.02 0.37 + 0.03 0533 490.02 
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S10, NA 37 .82 NA 
A1 503 4.30 + 0.14 20.4 + 0.2 4.05 + 0.07 
Feo0, 6.76 + 0.19 6.49 + 0.24 1.59 + 0.13 
MgO 1.03 + 0.11 0.36 + 0.28 1.86 + 0.11 
Ca0 125) = 058 3.04 + 0.13 27.9 + 0.5 
Na,0 0.423 + 0.008 0.344 + 0.006 0.115 + 0.003 
K50 1.56 + 0.21 3.77 + 0.19 1.68 + 0.12 
Tid, 0.25 + 0.04 1.06 + 0.06 0.19 + 0.04 
C org. NA NA 4.4 
As 29.8 + 0.4 87.1 + 0.7 8.03 + 0.25 
Ba 348 + 37 221 + 22 653 + 47 
Cl 628 + 30 469 + 22 180 + 4 
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Cu NA 49 NA 
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Mo 50 + 2 30a 8+1 
Ni NA 51 NA 
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Yb RE) oie 4.25 + 0.31 igp4 2 0M3 
Lu 0.38 + 0.03 0.43 + 0.03 Gets 0202 
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Si0, 48.44 39.02 57.74 
Al 503 aT 8.26 + 0.12 9.87 + 0.15 
Feo0, 10.5 + 0.2 2.81 + 0.18 2.38 + 0.17 
MgO 1.21 + 0.19 1.80 + 0.16 ND 
Cad Lees: i Tile 004 5.59 + 0.18 
Na,0 0.166 + 0.004 0.166 + 0.004 0.556 + 0.010 
K50 7.02 + 0.77 2.63 + 0.18 2.75 + 0.22 
Tid, 0.63 + 0.05 0.39 + 0.05 0.52 + 0.06 
C org. 1238 BA 4.1 
As 74.1 + 0.6 9.22 + 0.28 18s7 + 004 
Ba 419 + 27 209 + 33 283 + 33 
C1 354 + 18 238 + 23 400 + 23 
Co 38.0 + 0.3 Gig =: U2 leis 3 (ee 
Cr Tie CRE S 46.9 + 1.3 Te ies 
Cs 7.09 + 0.18 5.23 + 0.16 4.71 + 0.15 
Cu 143 43 69 
Ga 46.1 + 1.4 12.6 + 0.4 12.6 + 1.4 
Hf 00 Onl 1.4 + 0.1 es Gri 
Mn 103 + 2 428 + 6 175 +3 
Mo 167 + 9 18 +1 45 + 2 
Ni 391 69 105 
Pb 49 19 22 
Rb 186 + 6 91 + 4 Sorte 
Sb 9.69 + 0.10 1.14 + 0.04 5.75 + 0.08 
Sc 12.0 + 0.2 9.18 + 0.11 9.94 + 0.12 
Sr 476 + 30 533 + 41 93 + 25 
Ta 0.88 + 0.09 0.78 + 0.06 0.87 + 0.07 
Th 16.0 © 201 7.86 + 0.72 13.1 + 0.9 
U COPE: 9.32 + 0.20 15.2 + 0.2 
V 523 + 9 95 +5 450 + 10 
Zn 216 46 208 
La 55.1 + 0.3 17.4 + 0.6 21.3 + 0.5 
Ce 83.9 + 0.8 34.1 + 0.6 36.8 + 0.4 
Sm 13.8 + 0.5 3.21 + 0.36 5.40 + 0.30 
Eu 2.93 + 0.09 0.78 + 0.07 1.07 + 0.07 
Tb 1.70 + 0.08 ND ND 
Dy 9.80 + 0.14 2.20 + 0.09 3.60 + 0.08 
Yb Aik) & aye inte (et 1.93 + 0.19 
iA ND 0.19 + 0.02 0.27 + 0.02 
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eS SUES IN ieee AIX oh eee OLE Kove ee ee USE ee 
Si0, 52.90 NA 59.42 
Al 03 15.1 + 0.2 Ade & Ty Ge) ok ei 
Feo03 5.42 + 0.22 7.17 + 0.24 5.90 + 0.23 
MgO 1-51 + 0.23 0.75 + 0.29 0.76 + 0.20 
Ca0 0.62 + 0.09 320 + 0516 0.48 + 0.09 
Nap0 0.731 + 0.013 0.322 + 0.007 0.389 + 0.008 
K50 3.16 + 0.24 4.47 + 0.36 4.18 + 0.25 
Tid, 0.56 + 0.06 1.07 + 0.06 0.67 + 0.06 
C org: Dah NA 13 
As 40.1 + 0.5 86.7 + 0.8 35.4 + 0.7 
Ba 799 + 58 185 + 32 361 + 45 
Cl 760 + 33 473 + 22 477 + 23 
Co 33.2 + 0.3 16.7 + 0.2 Son ts ee 
Cr Mee 2 77.5 + 1.5 102 + 2 
Cs 8.56 + 0.19 5.19 + 0.18 Ie sip te One 
Cu 58 NA 53 
Ga 26.4 + 1.4 aoe ed 23.9 + 1.0 
Hf 3.6 + 0.1 9.5 + 0.2 We Ea 
Mn 96.6 + 2.6 75.4 + 2.1 122 + 3 
Mo 65 +2 2 ol ive] 
Ni 254 NA 107 
Pb 27 NA 77 
Rb 162 + 5 81+ 4 178 + 5 
Sb 10.4 + 0.2 2.67 + 0.07 6.97 + 0.13 
Sc Nae SS 2 15.4 + 0.2 ash ae (0.2 
Sr 183 + 26 107 + 23 112 + 24 
Ta 0.82 + 0.09 2.24 + 0.09 1.05 + 0.13 
Th 13.4)% 180 32.3 + 1.4 14.8 + 0.7 
U el se (06 21.8 + 0.3 5.55 + 0.15 
V 963 + 13 109 + 5 429 +7 
Zn 388 NA 269 
La 28.8 + 0.7 65.2 + 3.6 34.9 + 0.3 
Ce 53.5 + 0.9 129 +1 63.6 + 0.6 
Sm 4.23 + 0.30 14.3 + 0.4 4.38 + 0.29 
Eu 0.78 + 0.70 2.53 + 0.09 0.98 + 0.70 
Tb ND 2.03 + 0.08 ND 
Dy 2.85 + 0.08 9.05 + 0.14 3.45 + 0.08 
Yb 1.98 + 0.21 3.49 + 0.28 2.51 + 0.25 
Lu 0.32 + 0.03 0.44 + 0.03 0.41 + 0.03 
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SL UU so oe eee eS ee ee alse? oo 
S105 NA NA 24.24 
Al 503 9.20 + 0.09 AS 25etH0.0G 1538-+.0..04 
Fe503 3.82 + 0.19 14 -30400.0 pers Oi ys) 
MgO 2.04 + 0.17 0.60 + 0.11 ND 
CaO 1853-43023 14.894 0-2 ASS ag UES 
Nao0 0.170 + 0.004 0.342 + 0.006 0.121 + 0.003 
K50 3.39 + 0.24 1.42 + 0.26 0.14 + 0.06 
Ti0, Os29et 0.03 0.29 + 0.04 ND 
Ceaond. 30 siesta! 8.0 
As LOGO) E023 ee) ae WAS 1 O0mts Oni20 
Ba S33etes/ 336 + 40 ble 14 
Cl apn Cag 61251228 (elt So (als 
Co 1525 40.1 14.4 + 0.1 3401 + 0.3 
Cr. 54.4 + 1.4 30.0 stale SA) a2 (a8: 
Cs 17200-40012 Walks a Wigs 0.86 + 0.10 
Cu NA NA 52 
Ga 14.93 0.8 1059°+91.0 455 +1023 
Hf sie) ce 10k Cet Uel Oesec Oa 
Mn 410 + 6 Cliches 5OZ at 10 
Mo 60Rte2 40 +2 23401 
Ni NA NA 65 
Pb NA NA 11 
Rb 125 + 4 47 +9 1 Sars 
Sb 2o235 te 0200 EG MO ne 5100 (estas oe DOS) 
Se 1001 6.905402 12 (farseldy Ge OE) 
ST, Shey) wae Se) 29 tao) 17954336 
Ta 0.49 + 0.10 053340207 ND 
Th ea) se Wek 14.4 + 0.8 229 ete Ueo4 
U 123 ceraUec 37.8 + 0.4 GL ozet 020 
V ND ND 170 + 4 
Zn NA NA ou 
La (dary ie Weds T0079 C22 te0 6c 
Ce SPigi/ os Wat! SZere) USS) sae (RS 
Sm 4.30 + 0.35 SST (Se) arin Ge Os) 
Eu 0.83 + 0.08 4.07 + 0.14 Oe SE) ae WSO) 
Tb 0.25 + 0.06 Zed lets 0.08 ND 
Dy 3610+ 0.07 Wap he Ie Le ete0.08 
Yb 1.86 + 0.18 Shih ap (ashe Nesresy Gat a 
Lu 0.27 + 0.02 0245) 470-03 Q.14>+ 0502 
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Appendix 2 


Sample ID. 
EX, 10 Bx and 26 Ex 
EX 
Ex and 7 UBx 
Bx, thehxo and 24 hx 
Ex, 24 Ex and 34 Ex 
EX and 27 Ex 
Ex and 17 Ex 
Bx 2eo ke and oo Lx 
EX. 6 Hx, 20 Bx and 36 EX 
Bx coe Sc Eeeeands.oo08 hx 
Ex and 33 Ex 
Ex and 40 Ex 
Exe and 67 Ex 
EX 
Ex 
Exee CoE X, oO) ex rance ee kx 
Ex 
Ex 
Ex 
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Sample Locations. 


Well Location. 
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